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SECTION  I 


INTRODUCTION 


The  overall  function  of  the  Multi-Rate  Digital  Autopilot  is  to  provide 
stability  augnentation  far  the  pitch,  yaw,  and  roll  axis  of  an  aerodynamic 
or  space  vehicle  and  improve  dynamic  response.  The  generic  multi-rate 
loop  structure  is  similar  for  most  applications.  The  selection  of  the 
sanpling  and  iteration  rates  is  dictated  by  the  critical  performance  spec¬ 
ifications,  inherent  system  characteristics,  and  design  technique 
employed  (i.e.,  direct  digital  design  or  emulation).  Sanpling  here  refers 
to  the  rates  (i.e.,  samples  per  second)  at  which  signals  are  coupled  to 
and  from  the  digital  elements  within  the  autopilot.  The  iteration  rates 
are  the  rates  at  which  the  digital  filters  or  algorithms  are  updated. 

Both  these  rates  are  ocrmonly  referred  to  as  the  sample  rates  for  the 
autopilot. 

The  general  multi-rate  autopilot  is  a  hybrid  system  that  contains  both 
analog  and  discrete  elements.  The  analog  elements  include  the  airframe 
dynamics,  actuators,  sensors,  and  anti-aliasing  prefilters.  Discrete  ele¬ 
ments  include  those  autopilot  functions  performed  by  the  microprocessor 
(e.g.,  digital  filtering  and  gain  scheduling). 

The  sample  rates  in  a  missile  autopilot  can  typically  range  frcm  10  to 
600  Hz.  A  given  missile  autopilot  might  erploy  rates  of  10  ,  20  ,  40  ,  80, 
and  240  HZ.  These  widely  separated  rates  can  significantly  affect  the 
autopilot  performance  and  stability.  These  effects  are  further  accen¬ 
tuated  by  the  fact  that  missile  autopilots  are  typically  designed  using 
conventional  analog  filters  and  then  converted  to  a  digital  structure 


using  a  discretization  technique  (e.g. ,  Tustin  transform).  This  approxi¬ 
mating  procedure  is  called  emulation.  Since  the  emulated  digital  autopi¬ 
lot  only  approximates  the  analog  autopilot,  emulation  by  itself  can 
significantly  affect  performance  and  stability. 

Before  the  effects  of  multi-rate  'sampling  and  emulation  can  be 
assessed,  the  multi-rate  loop  structure  of  the  digital  autopilot  must  be 
modeled.  It  is  well  known  that  conventional  analog  or  single-rate  tech¬ 
niques  cannot  fully  model  the  multi-rate  autopilot  as  it  is  actually 
implemented  within  the  missile  microprocessor.  A  comprehensive  multi-rate 
methodology  is  needed  that  can  provide  a  high  fidelity  model.  This  model 
must  include  the  salient  features  of  the  multi-rate  sampling  scheme,  data 
transfer  and  storage  within  the  microprocessor,  and  actual  digital  filters 
implemented  in  the  software. 

This  report  presents  a  transfer  function  based  methodology  for  accu¬ 
rately  modeling  and  analyzing  the  multi-rate  digital  autopilot.  This  is  a 
'how  to'  report  that  offers  practiced  techniques  and  procedures  that  can 
readily  be  applied  by  the  control  engineer.  The  multi-rate  methodology  is 
applied  directly  to  a  typical  missile  autopilot.  However,  this  methodol¬ 
ogy  is  applicable  to  any  general  multi-rate  control  system. 

The  missile  autopilot  discussed  in  the  report  utilizes  angular  rate 
and  linear  acceleration  feedback  in  the  pitch  and  yaw  axis.  The  pitch  and 
yaw  channels  command  normal  and  lateral  missile  acceleration,  respec¬ 
tively.  Both  roll  and  roll  angular  rate  are  employed  to  maintain  a  rela¬ 
tively  constant  roll  axis  orientation. 

The  multi-rate  modeling  methodology  formulates  z-plane  equations  that 
are  well  suited  to  multi-rate  linear  analysis.  Time-domain  performance 


can  be  assessed  by  applying  a  sanpled  step  input  to  the  closed-loop  dif¬ 
ference  equation.  Effective  time  constant,  rise  time,  peak  time,  settling 
time,  and  percent  of  overshoot  can  be  determined  from  the  sanpled 
response.  In  the  frequency  domain,  bandwidth,  steady  state  gain,  phase 
margin,  and  gain  margin  are  typically  calculated  front  the  multi-rate  model 
equations. 

The  report  is  organized  into  five  major  sections  and  three  appendices. 
Sections  II  and  III  review  classical  sanpled  data  techniques  for  single¬ 
rate  and  multi-rate  control  systems.  This  review  starts  with  the  z- 
transform  and  progresses  to  multi- rate  systems.  Qrphasis  is  placed  on  the 
practical  application  of  multi-rate  analysis  techniques.  Detailed  mathe¬ 
matical  derivations  are  emitted  in  these  sections.  Appropriate  references 
are  cited  that  cover  the  detailed  analytical  theory. 

Section  IV  surrmarizes  the  specific  multi-rate  autopilot  that  is  used 
throughout  the  report  to  illustrate  the  multi-rate  modeling  methodology. 

In  Section  V,  the  multi-rate  methodology  is  applied  to  the  pitch  channel 
of  the  autopilot.  Both  open- loop  and  closed- loop  multi-rate  equations  are 
generated.  Section  VI  presents  three  computer  programs  that,  together, 
can  be  used  to  numerically  evaluate  the  multi- rate  model  equations.  The 
Appendices  expand  the  material  in  the  report  and  present  illustrative 
examples  to  help  clarify  the  multi-rate  methodology. 


m 
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SECTION  II 

REVIEW  OF  FUNDAMENTAL  Z-TRANSPORM  THEORY 

A.  INTRODUCTION 

The  fundamental  principles  of  the  z- transform  are  presented  in  this 
section  (References  1  through  10).  Practical  z- transform  theory  is  based 
on  the  assumption  that  the  actual  sampling  operation  can  be  modeled  as  the 
amplitude  modulation  of  an  inpulse  train.  This  central  concept  greatly 
simplifies  the  analysis  and  synthesis  of  sampled  data,  discrete,  or  digi¬ 
tal  systems.  This  view  of  the  sampling  process  is  valid  for  most  prac¬ 
ticed  systems  and  use  of  this  theory  is  "normally  considered  exact. 

The  multi-rate  digital  autopilot  contains  both  continuous  and  discrete 
elements.  The  z- transform  provides  a  unified  analysis  and  synthesis  tech¬ 
nique  far  these  hybrid  systems.  For  a  sampled  continuous  element,  the  z- 
transfcon  can  be  considered  as  the  Laplace  transform  of  an  impulse 
sequence  (inpulse  train)  where  the  area  or  strength  of  each  individual 
impulse  equals  the  value  of  the  continuous  time  function  at  a  particular 
discrete  sampling  instant.  An  alternate  viewpoint  is  to  consider  the 
exponent  in  the  z~n  delay  operator  as  an  ordering  variable  for  a  number 
sequence  (or  a  sequence  of  discrete  signal  values)  where  the  coefficient 
far  the  z“n  term  equals  the  value  of  the  number  sequence  (or  the  discrete 
signal)  at  the  nth  discrete  time  delay.  This  second  viewpoint  is 
exemplified  by  the  direct  modeling  of  the  time-dcmain  difference  equation 
in  the  z-plane  ( see  Section  II-F) . 

In  practice,  to  obtain  a  discrete  z-plane  model  of  a  hybrid  control 
system,  the  continuous  time  functions  are  first  expressed  in  the  s-plane 
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and  then  transformed  to  the  z-plane  using  standard  techniques  such  as  par¬ 
tial  fraction  expansion  coupled  with  table  look  up.  A  second  z- transform 
method  employs  the  inversion  integral  and  contour  integration  (Reference 
4) .  The  partial  fraction  expansion  approach  is  mechanized  in  the  DISCRET 
computer  program  ( Reference  11).  Discrete  functions  (e.g. ,  digital 
algorithms  in  a  microprocessor )  are  first  modeled  with  a  difference 
equation  and  directly  converted  to  the  z-plane  by  substituting  the 
appropriate  z“n  delay  operator  for  each  term  in  the  difference  equation. 
During  the  design  phase,  it  is  the  z-plane  filter  that  is  first  obtained. 
This  filter  is  then  converted  to  a  difference  equation  using  the  z“n  delay 
operator  and  subsequently  implemented  on  a  digital  microprocessor. 

No  distinction  is  necessary  between  the  z- transform  function  derived 
from  a  sanpled  continuous  element  and  the  z-transform  function  that  models 
an  inherently  discrete  element.  Once  discretized,  all  elements  of  a 
hybrid  system  can  be  treated  using  ccnrnon  analysis  and  design  techniques. 
However,  the  standard  z-transform  models  continuous  functions  at  their 
sanple  points  and  does  not  provide  inter- sairple  information.  It  may  be 
necessary  to  investigate  the  inter- sanple  behavior  of  discretized  con¬ 
tinuous  elements  using  such  techniques  as  the  advanced  (or  delayed)  z- 
transfarm  or  the  “continuous  frequency  response  of  a  discretely  excited 
system"  (Reference  12) . 

B.  DEFINITION  OF  THE  z-TRANSFORM 


The  definition  of  the  z-transform  stems  fran  the  infinite  surmation 


C?(t)  =  c(0)  +  c(T)  6 ( t-T)  +  c(2T)  S(t— 2T)  +  ...  +  c(  )  S(t-  ) 

-C  c(kT)  S(t-kT)  ,  k  =  0,  1,  2,  ...  (1) 


(2) 


An  alternate  expression  for  Equation  1  is 


<^(t)  =  c(t) 


OP 

y  &(t-y 


=  c(t)  6  T(t) 

The  inpulse  train  6  T(t)  represents  a  series  of  inpulses  of  unit  strength 
or  area  equally  spaced  in  time  and  extending  fran  zero  to  plus  infinity. 

In  Equations  1  and  2,  the  sampled  signed  (^(t)  is  represented  by 
inpulses  whose  area  or  strength  equals  the  magnitude  of  the  continuous 
signal  c(t)  at  the  sanpling  instance  t=kT.  This  view  of  the  sampling  pro¬ 
cess  as  the  amplitude  modulation  of  an  irrpulse  train  6  T(t)  by  the  con¬ 
tinuous  signal  c(t)  (Figure  1)  farms  the  mathematical  basis  fear  the  z- 
transfarm.  Such  a  viewpoint  is  justified  if  the  actual  time  during  Which 
the  sanpler  is  closed  is  short  compared  to  the  time  constants  in  the 
system  under  investigation.  It  is  shown  in  Reference  5  that  for  a  system 
with  a  single  time  constant  tau  =  l/a,  the  error  using  impulse  modulation 
is  less  than  5  percent  for  a  sampler  pulse  width  h  Which  is  less  than  or 
equal  to  one-tenth  of  the  time  constant  ( i.e. ,  h/tau  <  1/10) . 


In  general,  if  the  Laplace  transform  of  the  continuous  time  function  c(t) 
is  a  rational  algebraic  function,  a  closed  form  can  be  found  for  the  infi¬ 
nite  series  representation  of  (^(s).  The  simple  change  of  variable 


Amplitude  Modulation  of  Impulse  Train 
by  Continuous  Signal 


(4) 


z  =  e8*1,  results  in  the  one-sided  z- transform. 

go 

<*.,  =  ^  c(kT)z~k 
k=0 


For  the  two-sided  z-transform,  the  lower  suimation  limit  becomes  minus 
infinity  and  c(t)  is  defined  for  negative  time.  Most  practiced  applica¬ 
tions  utilize  the  cne- sided  z-transform. 


C.  TABLE  OF  z- TRANSFORMS 


Utilizing  the  fundamental  definition  of  the  z- transform  in  Equation  4, 
Reference  1  generates  a  table  of  z- transforms  for  elementary  continuous 
time  functions.  This  table  is  repeated  here  as  Table  1.  The  column 
labeled  F(z,  A  )  is  for  the  advanced  z-transform  with  0  <  A  <1.  More 
extensive  z- transform  tables  can  be  found  in  References  2  through  8. 

The  z- transform  variable  in  Table  1  is  defined  as  z  =  e57,  where  t  is 
the  sampling  interval  in  seconds.  For  other  sanpling  intervals  such  as 
T/3,  we  sinply  replace  T  with  T/3  in  the  z-transform  expression.  For 
example,  consider  the  continuous  time  function  g(t)  =  te~at.  The  Laplace 
transform  for  this  function  is  l/(s  +  a)  2  and  the  s-plane  sampled  function 
is  given  by 

dT/3(s)  „  c - i - .]  (5) 

(s  +  a)2 


The  z- transform  far  Equation  5  can  be  obtained  directly  fran  Table  1. 


GT/3(Z,  .  (T/3)ze-*T/3 
(z  —  e”aT/2)2 


(6) 


where 


f(t)  F(b)  F(*)  F(z,  A)  =  Z[f(t  +  AT)] 


coa  bT)  +  (o  fin  bT)  ■  «  -  So  cob  W*  +  a 


dr/3(z3)  =  (T/3 ) z3e~aT/3 

{Z3  -  e“aT/3)2 


where 

Z3  =  esT/3  (9 

The  T/3  superscript  cn  the  function  GT/3(s)  and  GT/3(z)  designates  the 
sanpling  interval  in  seconds.  The  use .of  the  superscript  notation 
becomes  necessary  in  multi-rate  systems  to  separate  and  specify  the 
sanpling  interval  for  each  functioned  element.  The  multi-rate  notation 
is  covered  in  Section  III. 

D.  PARTIAL  FRACTION  EXPANSION  AND  TABLE  LOCKUP 

Given  a  sanple  continuous  time  function  expressed  as  a  sampled  s- 


plane  transfer  function  (e.g. ,  GT(s)),  partial  fraction  expansion  and 
table  lockup  can  be  used  to  obtain  the  z-plane  transfer  function.  The 
sampled  system  in  Figure  2  will  be  used  to  illustrate  this  approach. 


The  output  equation  for  Figure  2  is  given  by 


d?  =  (gm)trt 


-  1  -  e-sT 


G(s)  = _ a 


s  +  a 


so  that 


(GM)t  =  r(l  -  e-sTaiT  = 


s(s  +  a) 


<* -«-*■>  '4  -  rhr>T 


M(s)  in  Equation  11  is  the  s-plane  transfer  function  for  a  zero-order-hold 
(ZOH).  In  Equation  12,  the  a/s(s  +  a)  term  is  separated  into  two  terms 
using  partial  fraction  expansion  (Reference  1) .  Since  z  =  esT, 

1  “  e_sT  =  1  -  z-1  *  (z  -  l)/z  and  this  z-plane  term  can  be  separated  from 
the  sanpling  operation  (see  Section  III-D) .  Therefore,  from  Table  1  we 


(GM)T  _  z  ~  1  ( — 5 —  -  - 5 - )  _ 

z  z  -  1  z  _  e-aT 


1  -  e~aT 
z  -  e-aT 


E.  THE  Zr-TRANSTORM  INVERSION  INTEGRAL 


The  inversion  integral  is  a  closed  form  technique  for  finding  the 


inverse  z- transform  (Equation  14) . 


c(kT) 


^c'T(z)  dz 


Equation  14  is  based  on  the  Laurent  series  expansion  of  F(z)  =  zk“^dr(z) 
about  z  =  0.  To  develop  Equation  14  we  first  expand  Equation  4,  the  fun¬ 
damental  definition  of  the  z- transform,  * 


cfT(z)  =  c(0)  +  c(T)z_1  +  c(2T)z“2  +  ...  +  c(W)z-k  + 


•  •  • 


(15) 


If  we  new  multiply  Equation  15  by  z*4-*,  we  get 


F*(z)  =  ^-IdT(z)  =  c(0)z*“l  +  c(T)z^“2  +  ...  +  c(l<T)z“l  +  ...  (16) 

The  desired  output  c(kT)  in  the  Laurent  series  expansion  (Equation  16)  is 
defined  as  the  residue  of  the  function  F^z).  To  obtain  Equation  14  from 
Equation  16,  we  take  the  closed  contour  integral  of  Equation  16  and  apply 
the  Cauchy  Theorem.  This  theorem  states  that  if  a  carpi  ex  integral  is 
defined  by 


1  2?^  dz  (17) 

2  it  j 

Cl 

and  the  integral  is  taken  around  a  closed  contour  q  which  encloses  the 
origin  of  the  z-plane ,  then  the  value  of  the  integral  reduces  to 

z*dz  =  2  tt  j,  for  k  =  -1  (18) 

*  0  ,  for  k  ^  -1 

Continuing  with  the  derivation,  the  closed  contour  integral  of  Equation  16 
is  given  by 


1 

2  it  j 


Cl 


z^^-d^z)  dz  =  1 

2  IT  j 

=  c(kT) 


Cl 


C(0)z^”^  +  ...  +  c(kT)z"l  +  ...dz 

(19) 


The  right  hand  side  of  Equation  19  reduces  to  c(kT)  since  only  the 
Z“1  term  the  integration  is  non-zero.  That  is, 

c(kT)  zT1  dz  =  yi-j  c(kT)  2  tt  j  =  c(kT)  (20) 

q 

The  left  hand  side  of  Equation  19  is  the  inversion  integral  given  in 
Equation  14.  Applying  a  second  Cauchy  Theorem,  this  integral  reduces  to 


the  evaluation  of  the  residues  of  the  poles  associated  with  the  integrand 
zk-lcfT(z) .  Equation  19  can  then  be  expressed  as 


c(KT)  *  ^ residues  of  Qf  2k-ldr(z) 


Consider  the  following  example  for  Equation  21.  Let 


Inserting  Equation  22  into  Equation  21  produces 


c(kT)  =  residues  zk“l  _ ,8z _ 

(z  -  .1)(  z  -  .« 


.8zk 
z  -  .1 


z  =  .9 


z  =  .1 ,  .9 


z  =  .1 


(•9)k  -  (.1)* 


F.  SAMPLED  OR  DISCRETE  TIME  RESPONSE 


The  discrete  time  response  for  a  z-plane  transfer  function  can  be 
calculated  using  the  z- transform  inversion  integral  (Equation  21).  The 
use  of  this  suimation  is  not  always  practical  or  convenient  since  it 
requires  the  calculation  of  residues.  Appendix  C  in  Reference  1 
demonstrates  various  techniques  for  obtaining  the  discrete  time  response 
This  subsection  presents  two  of  these  techniques,  the  recursion  or  dif¬ 
ference  equation  and  continued  fraction  expansion.  The  single-rate 
sampled  system  in  Figure  3  will  be  used  to  illustrate  these  two  tech- 


Figure  3.  Single-Rate  Sampled  System 


In  the  continued  fraction  expansion  approach,  the  output  sampled 
variable  d?(z)  is  expressed  as  a  product  of  the  z-plane  transfer  function 
for  the  system  and  the  z-plane  expression  for  the  sampled  input. 

CF(Z)  =  GT(z)RT(z)  (24) 

The  nuneratar  of  this  z-plane  product  Gt(z)Rt(z)  is  divided  by  the  denani- 
natar  and  the  negative  powers  of  z-n  in  the  quotient  are  interpreted  as 
impulse  functions  shifted  in  time  [i.e. ,  z-n  =  e-snT  =  6  (t-nT)] .  Suppose 

_  «75z 

W  "  z*  -  ,75z  +  .125  ,  Rt(z)  =  1  (25) 

Then 


C?(z)  =  (jT(z)rT(z) 


.75z 

z2  -  ,75z  +  .125 


=  .75zrl  +  .5625z'2+  .328125z"3  +  .17578125Z-4  +  ...  (26) 


Frcm  the  negative  exponents  in  Equation  26,  the  sampled  output  is  written 
directly  as 


Option  66  in  the  TOTAL  computer  program  (Reference  13)  divides  polynomials 
and  can  be  used  to  implement  the  continued  fraction  expansion  technique. 


The  recursion  or  difference  equation  approach  provides  a  closed  form 
expression  for  the  discrete  time  response.  We  first  write  the  pulse 
transfer  function  relating  the  sampled  output  to  the  sampled  input.  For 
exanple,  rewriting  Equation  26  gives 

cF(z)  =  Gt(z)  =  ,75z 

Rt(z)  z2  -  .75z  +  .125  (28> 


The  pulse  transfer  function  is  next  converted  to  negative  exponent  form  by 
multiplying  the  numerator  and  denominator  by  z-n  where  n  is  the  highest 
power  of  z  that  appears  (z2  in  Equation  28).  Equation  28  then  becomes 


C?(z) 

Rt(z) 


_ .75z~x 

1  -  «75z_1  + 


.125z“ 


(29) 


Cross  multiplying  the  numerators  and  denominators  in  Equation  29  produces 

C?(z)  Cl  -  .75z_1  +  .125z-2]  =  Rt(z)  C.75Z"1]  (30) 

To  write  the  difference  equation,  the  negative  powers  of  z  are  interpreted 
as  delayed  impulse  functions.  That  is, 


Z-1  _  e-sT  -  5  (t-T) 

2-2  -  e-2sT  -  6  (t-2T) 


z-n  =  e-nsT - 5  (t-nT) 


(31) 
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c(kT)  =  .75c  [(k-l)T]  -  .125c  [(k-2)T]  +  .75r  [(k-l)T] 


(32) 


Equation  32  indicates  that  the  present  value  of  the  sartpled  output  c(kT) 
is  a  function  of  the  past  two  values  of  the  output  and  the  past  sampled 
value  of  the  input.  This  closed  form  equation  is  calculated  by  Option  35 
in  TOTAL  (Reference  13). 


Since  the  outer  sample  rate  operates  through  the  inner  rate  (see  subsec¬ 
tion  C  of  this  section) ,  Equation  34  reduces  to 

dT/N  =  GT/NRT  (35) 

The  superscript  notation  in  Equation  35  designates  the  time  interval 
between  each  sampling  operation.  For  example,  gt/n  indicates  that  the  func¬ 
tion  G  is  sampled  at  N/T  sanples  per  second  ano  RT  denotes  that  the  signal 
R  is  saitpled  at  l/T  sanples  per  second.  Also,  if  GT/N  represents  a  digi¬ 
tal  filter  implemented  in  a  microprocessor,  then  T/N  is  the  update  inter¬ 
val  of  this  filter  in  the  software.  These  two  interpretations  for 
gT/N  will  be  discussed  later.  The  traditional  notation  for  RT  and  G>r/N  is 
R(z)  and  Giz^) ,  where  z  =  esT  and  z^  =  esT/N.  However,  the  superscript  nota¬ 
tion  is  used  in  the  report  because  it  permits  an  explicit  statement  of  the 
sairpling  rates  with  minimal  notational  complexity. 

The  general  notation  GT/^  implies  either  GT/N(s)  or  GT/N(zn) ,  they  are 
equivalent  in  the  discrete  time  domain.  GT/N(s)  represents  the  Laplace 
transform  associated  with  the  sampled  continuous  time  function  gCkT/N) . 

G(s)  is  the  Laplace  transform  of  the  continuous  function  g(t) .  GT/N(s)  is 
given  by  the  standard  definition  of  the  z- transform  (see  Section  II) . 

O O 

gT/N(s)  =  H  g(kT/N)e_ksT/N  (36) 

k=0 

The  equivalent  <jI’/N(zn)  function  is  the  z- transform  expression  obtained 


from  with  =  esT/N.  That  is. 

oo 

GT/N(Zn)  =  GfT/N(s) 

y  g(hT/N)z~k 

z  =  esT/N 

k=0 

(37) 

This  simple  change  of  variable  {z^  -  es'I'/N)  converts  GT/N(s),  a  nonalgebraic 


function  in  s,  to  a  rational  function  in  z. 

A  second  interpretation  for  the  GT/N( 7^)  function  is  also  possible. 
G^/^Zn)  could  represent  the  z-plane  model  of  an  inherent  discrete  function 
(e.g. ,  digital  filter  implemented  in  a  microprocessor).  In  this  case, 
is  directly  associated  with  a  discrete  time  domain  difference 
equation  and  is  not  associated  with  a  sampled  continuous  time  function 
g(kT/N)  as  is  true  in  Equation  37. 

To  help  clarify  the  equivalence  between  GT/N(s)  and  GT/N(zn)  and  to 
illustrate  the  two  possible  interpretations  for  the  z-plane  function 
G^/N(zn)#  consider  the  single  rate  system  in  Figure  6. 


Figure  6.  Single-Rate  Sampled  System 


The  output  equation  far  Figure  6  is 


C?  =  Gt(s)Rt 


1  -  e-sT  1 
s  s  +  1 


Since  z  =  esT, 


C/T(s)  =  (1  -  z-1)  (-,  T 


sTs  +  f] 


=  ^  a  “  s"TT* 


Fran  the  table  of  z- transform  in  Section  II,  the  GT(z)  function  for 
Equation  41  is  given  by 


GT(z)  =  (1  -  z-1)  ( 


z  -  1  z  -  e~T 


= 


(1  -  e~T)z 
{z  -  1)  (z  -  e“T) 


1  -  e-T 

z  —  e~ ^ 


$3  * GT(rf 


(1  -  e~T)z~^ 
1  -  e  -Tz-1 


The  difference  equation  for  Equation  43  is  obtained  by  cross  multiplying 
the  numerators  and  denominators  and  interpreting  z“l  as  a  one  sample 
period  delay. 


C(kT)  *  e-TC  [(k-l)T3  +  (1  -  e“T)R  [(k-l)T] 


Equation  44  is  the  difference  equation  for  the  z-plane  function  GT(z)  in 
Equation  -»J  Which  is  equivalent  to  the  sampled  s-plane  function  GT(s)  in 
Equation  40. 


A  second  interpretation  for  the  general  GT(z)  function  can  be  obtained 


from  the  difference  equation 

C(kT)  =  .9CC (k-l)T]  +  R(kT)  +  .5R  [(k-l)T]  (45) 


Suppose  Equation  45  represents  the  code  in  a  microprocessor  vhich  imple¬ 
ments  the  digital  filter  given  by 


G?(z) 


z  +  .5  _  1  +  «5z~l 
z  -  .9  ~  1  -  .9z_1 


(46) 


In  this  case,  the  GT(z)  function  in  Equation  46  models  the  discrete  time 
domain  operations  in  a  digital  microprocessor;  it  is  not  directly  i elated 
to  a  sanpled  continuous  function  (as  is  true  far  Equation  43). 


In  multi-rate  system  design  and  analysis  no  distinction  is  required 
between  the  z- transform  function  GT/N(z)  derived  from  a  sampled  continuous 
function  (e.g.,  )  ard  the  z- transform  function  that  models  the 

digital  operations  in  a  microprocessor.  The  z- transform  (and  its  exten¬ 
sions  the  w-  or  w‘  -transform)  provides  a  unified  analysis  and  synthesis 
technique  far  systems  containing  both  continuous  and  discrete  elements 
( i.e. ,  hybrid  systems).  This  ability  to  model  continuous  and  discrete 
elements  in  a  common  discrete-domain  is  perhaps  the  most  fundamentally  use¬ 
ful  property  of  the  z- ,  w-,  and  w' -transforms. 


C.  EQUIVALENT  HIGH- RATE  MODEL  FOR  MULTI-RATE  SYSTEMS 

The  general  multi-rate  transform  product  illustrated  by  Equation  35 
(repeated  below  as  Equation  47)  is  addressed  in  this  subsection. 


C^/^Zn)  =  <jr/N(zn)RT(z) 


(47) 


Equation  47  cam  be  placed  in  closed  form  by  establishing  a  cannon  high- 
rate  definition  for  the  z- transform  variable.  That  is, 

dr/N(zn)  =  tfr/N^R^z)  ,  GT/N(Zn)RT(zN) 

“  /  4/^ 


where 


s  e8*11,  Zn  =  esT/N 


To  evaluate  Equation  48,  we  first  calculate  the  T/N  and  T  z- transforms  for 
GT/N(zn)  and  Rt(z),  respectively.  The  trams form  variable  z  in  RT(z)  is 
next  replaced  with  z  =  z^  .  The  transform  product  GT/N(zn)RT(z^) 
cam  then  be  obtained  since  appears  in  both  functions. 

Consider  an  example.  Let 


*■>  ■; 


a(3>  -  ?Ta 


^ n)  -  <-§>t 


Referring  to  the  z-transform  table  in  Section  II, 


For  N  =  2  we  have 


dr/2<^> ' 


(Z2  -  e"aT/2)  (z2  -  1) 


-  e  -  Z2  +  e  * 


(53 


D.  AN  IMPORTANT  MULTI-RATE  TRANSFORM  RULE 

It  is  well-established  (References  1-12)  that  for  the  single-rate 
system  in  Figure  6,  the  output  equation  is  given  by 

C  =  CRT  (54 

Then 

C?  =  (GRT)T  =  GTrT  (55 

and 

d?  =  Gt  (56 

RT 

or 

C?(z)  _  Gt(z),  z  =  e37  (57 

RT(z) 

In  Equation  57,  z)  is  defined  as  the  pulse  transfer  function  relating 
the  sampled  output  C^(z)  to  the  sanpled  input  RT(z).  We  observe  that 
the  RT  discrete  function  in  Equation  55  is  simply  separated  fran  the 
resampling  operation,  as  was  done  in  Equation  35.  These  two  examples 
illustrate  a  more  general  rule. 

The  outer  sampling  operator  (  .  )T/N  operates  through 
an  inner  sampling  operator  if  the  ratio  of  the  inner 
sampling  period  to  the  outer  sampling  period  is  an 
integer. 


If  a  multi-rate  transform  product  is  resampled,  discrete 
functions  within  the  product  can  be  separated  frcm  the 
resanpling  operation  if  the  resampling  rate  is  equal  to 
or  greater  than  the  sanple  rates  of  the  individual  dis¬ 
crete  functions. 

This  general  multi-rate  transform  rule  is  illustrated  in  Equation  58  (N 
and  M  are  positive  integers). 

UG1G2G3  ...  Gi)  (rTrT/2  ...  rT/NrT/W-1  ...  rT/WM)]  T/N  = 

C (G1G2G3  ...  GiMRT/fl  ...  rJ/^M)]  T/N  I^rT/2  ...  rT/N  (58) 

The  C(G1G2G3  ...  Gi)(R^/Nt1  ...  R^AjfM)]T/N  transform  operation  represents 
high-to-lcw  rate  transform  oonversion  since  the  resample  rate,  N/T  samples 
per  second,  is  lower  than  the  inner  discrete  transform  rates.  The  func¬ 
tions  G]G2G3  ...  Gj[  are  e-plane  continuous  functions  and  ... 

are  discrete  functions.  (A  numerical  example  of  a  high- to- low  irate  trans¬ 
form  oonversion  is  presented  in  Appendix  B. ) 

An  explicit  example  that  illustrates  the  multi-rate  transform  rule  in 
Equation  58  is  given  in  Equation  59. 

CG^(s)G2(s)G^(z)G^/^(Z2)g^/2(z3)G^/^(z4)]  = 


CGi(s)G2(s)OT/3(Z3)GT/4(z4)]  T/2  gT(z)gT/2(z2) 


(59) 


Equation  101  fran  Section  V  also  applies  the  multi-rate  transform  rule  in 
Equation  58.  This  example  is  repeated  below. 


ikzL  =  [SUl  SCT(s)m2(s)gJ/2(22)gJ/2(z2)gJ(z)]T 
A2c(z) 

«  [Mil  ACT(s)M2(s)gT/2(Z2)GT/2(Z2)]T  gT(z) 

=  [(Mil  ACTCsjM^s))^2  gT/2(z2)GT/2]T  gt(z)  (61) 

The  numerical  evaluation  of  the  high- to- low  rate  transform  conversion  in 
Equations  59  and  61  requires  a  phantom  sampler  and  the  calculation  of 
residues.  In  Equation  61,  the  phantom  sampler  is  used  to  sample  the  s- 
plane  functions  at  a  T/2  interval  prior  to  calculating  the  high- to- low 
rate  conversion.  The  phantom  sampler  principle  and  the  high- to- lew  rate 
transform  conversion  are  covered  in  the  next  subsections. 

E.  HIGH-TO-LOW  RATE  SAMPLING  AND  THE  PHANTOM  SAMPLER 

Consider  the  fundamental  fast- input/slow-output  multi- rate  sampling  con¬ 
figuration  in  Figure  7.  M  represents  the  transfer  function  of  a  data  hold 
and  G  a  continuous  system  in  the  s-plane.  In  general,  the  output  sampling 
interval  T^  is  greater  than  the  input  sampling  interval  T^j  however,  a 
direct  integer  ratio  between  output- to- input  sanpling  is  not  necessarily 
inplied.  That  is,  Tj  is  not  necessarily  equal  to  T^  times  an  integer 
(e-g. ,  T2  a  6T  ).  This  simple  multi-rate  configuration  represents  a  prao- 


tical  situation  which  can  be  analyzed  using  a  cannon  sanpling  period  T  such 
that 


T  _  _  T  „  M,N  =  integer  (62) 

M  l  N  2 


Figure  7.  Fundamental  Fast-  Input/ Slew-Output 
Multi-Rate  Sampling,  T2  >  Ti 


For  exanple,  if  the  input  is  sampled  at  20  Hz  and  the  output  at  13.33  Hz? 
=  0.050,  T2  =  0.075,  and  a  choice  of  T  =  0.150  produces 

|  -  TL  I  =  T2  T  «  0.150  (63) 

The  multi-rate  configuration  in  Figure  7  can  be  reduced  to  the  general 
farm  in  Figure  8. 


Figure  8.  Fast-Input/Slow-Output  Sanpling 
with  Ccrrmon  Sanpling  Period  T 


The  output  equation  for  Figure  8  is  given  by 
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(64) 


cT/N  =  Cgmrt/m] 


Computationally,  Equation  64  is  more  involved  than  a  general 
slew- input/ fast-output  system,  since  the  larger  T/N  sanpling  interval  does 
not  operate  through  the  snaller  T/M  sanpling  interval  (see  subsection  D 
of  this  section) .  Equation  64  can  be  evaluated  by  mathematically  adding  a 
phantom  sanpler  to  the  output,  Figure  9.  The  sarrple  rate  for  the  phantom 
sampler  must  be  an  integer  multiple  of  the  output  sample  rate.  This 
mathematical  trick  is  valid  since  the  slower  rate  output  sanpler  (N/T 
sanpler  per  second)  simply  rejects  all  the  unwanted  samples  from  the 
higher  rate  phantom  sanpler  (l/T*  sanpl.es  per  second).  With  the  addi¬ 
tional  phantom  sanpler  T*,  the  output  equation  far  Figure  9  becomes 

< ?'»  -  C(®)T*rT/M]T/N  (65) 


Figure  9.  Fast-Input/Slow-Output  Sanpling  with 
Phantom  Sanpler  T* 


The  problem  new  reduces  to  selecting  T*,  calculating  the  transform  pro- 
T  T /M 

duct  GM  *R  '  ,  and  calculating  the  high- to- lew  rate  transform.  As  pointed 
out,  T*  must  be  an  integer  suhmultiple  of  the  output  sanpling  interval  t/n. 
In  addition,  T*  must  have  an  integer  relationship  with  the  sanpling  interval 
T/m  to  allow  the  inner  transform  product  to  be  calculated.  For  the  most 
general  ca^e  where  m/n  in  Figure  9  is  not  an  integer,  an  obvious  choice  for 
the  phantom  sanpler  is  T*  =*  T/f*I.  However,  if  M/N  is  an  integer,  the  best 


choice  is  T*  =  T/M.  For  example,  if  M  =  3  and  N  =  2,  T*  =  t/6  and 
Equation  65  becomes 

CT/2(z2)  =  [(GM)T/6(z6)RT/3(z3)]  ^  (66) 

2 

Since  T/6  and  T/3  have  an  integer  relationship,  Z3  —  and  the  inner  trans¬ 
form  product  cam  be  calculated  using  z&  as  a  cannon  variable.  That  is, 

(GM)T/o(z6)RT/3(z3)  =  (GM)T/6  (z6)RT/6(z6)  =  Gj/6(z6)  (67) 


T /6  2 

where  the  R  '  (z^)  discrete  function  is  formed  by  simply  substituting  zg 

T/3 

for  every  z3  in  the  R  '  (z^)  discrete  function.  Making  this  substitution. 
Equation  66  reduces  to  the  general  form 

cT/2(z2)  =  [gJ’/6(Z6)]  t/2  (68) 

where 


*2 


a  esT/2 


(69) 


The  high- to- lew  rate  transform  conversion  in  Equation  68  can  be  calculated 
using  Equation  71  in  subsection  F.  Note  that  the  use  of  T*  =  T/MN  =  t/6 
has  allowed  the  formation  of  the  inner  transform  product  and  at  the  same 
time  provided  an  integer  ratio  for  the  outer- to- inner  sampling  intervals. 


F.  HIGH-TO-LCW  RATE  TRANSFORM  CONVERSION 


A  computer  program  7XC0NV  is  presented  in  Reference  11  that  calculates 
a  lew-rate  discrete  transform  frem  a  given  high-rate  discrete  transform. 
The  general  form  of  this  transform  conversion  is  given  by 

cTU)  -  CcT/N(zn)]T  .[GT^RT/H]T  (70) 


where 
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C  (z)  =  low-rate  z- plane  function,  z  =  e 

CT/N(zn)  =  high- rate  z-plane  function,  zn  = 

This  tr ans formation  arises  vhen  a  high-rate  system  element  in  a  multi-rate 
control  system  is  resampled  at  a  lower  rate  (Figure  10) .  The  high- to- low 
rate  conversion  in  Equation  70  can  be  evaluated  using  Equation  71  (see 
Reference  11)  far  a  derivation  of  Equation  71. 

zn  =  poles  C^/N(zn)  (71) 
*n 

The  TXCONV  computer  program  mechanizes  the  residue  calculations  in 
Equation  71.  TXCONV  is  also  available  as  options  144-149  in  the  TOTAL 
computer  program  (Reference  13) .  The  user  simply  supplies  the  high-rate 
z-plane  discrete  function  GT/N(zn)RT/N(zn)  and  options  144-149  calculate 
the  dT(z)  discrete  function.  Appendix  B  presents  an  example  in  which  the 
residue  calculations  in  Equation  71  are  illustrated. 


SECTION  IV 


MULTI-RATE  AUTOPILOT  BLOCK  DIAGRAM 

A.  OVERALL  BLOCK  DIAGRAM 

The  block  diagram  in  Figure  11  depicts  the  specific  multi-rate  autopi¬ 
lot  structure  used  throughout  the  report.  Angular  rates  (p,  q,  and  r)  and 
linear  acceleration  (Az  and  Ay)  frcm  body  mounted  sensors  provide  the 
feedback  signals  for  the  autopilot.  The  roll  angle  0  is  obtained  by 
integrating  the  output  frcm  the  roll  rate  gyro.  Angular  rate  feedback 
augnents  the  low  natural  aerodynamic  dartping  of  the  airframe. 

Acceleration  feedback  provides  the  required  steady  state  and  transient 
acceleration  response  to  ccrrtnand  inputs.  and  Ayc  represent  ccmnand 

inputs  to  the  pitch  and  yaw  autopilot  loops,  respectively.  Hie  digital 
update  rates  are  shewn  as  superscripts  with  each  digital  filter.  For 
example,  the  T/2  superscript  indicates  that  the  (ATF)T/2  digital  filter  is 
updated  every  T/2  seconds. 

(ATF)t/2  in  the  autopilot  is  a  digital  notch  filter  that  suppresses  the 
first  body  bending  mode  of  the  flexible  airframe.  (KTF)t  is  a  digital  lead 
filter  that  improves  the  overall  missile  acceleration  response  to  ccmnand 
inputs.  ACT  is  the  actuator  transfer  function,  and  A/C  is  the  missile 
airframe  transfer  functions.  Analog  filter  models  the  storage  of  data 
within  the  microprocessor.  M2  models  the  digital- to- analog  conversion  pro¬ 
cess.  Hie  remaining  digital  filters  provide  signal-conditioning 
( compensation)  to  improve  missile  response  and  stability. 

B.  PITCH/YAW  AXIS  AUTOPILOT 

The  autopilot  in  Figure  11  contains  identical  pitch  and  yaw  loops 


(ITF) 


t  Block  Diagram 


{ symmetrical  airframe  dynamics).  The  pitch  channel  ccrtnands  effective  ele¬ 
vator  deflections  (  6  e)  and  the  yaw  channel  rudder  deflections  (  6  r) .  In  the 
pitch  loops,  body  angular  rate  (q)  abcxrt  the  y-axis  is  first  sampled  every 
T/6  seconds,  averaged  over  three  samples,  resampled  at  a  T/2  interval,  and 
then  feedback  through  filters  (DTF)T/2  and  (ATF)T/2.  This  three  sample 
averaging  is  modeled  with  digital  filter  (JTF)'1'/®  (see  subsection  IV-E) . 
Linear  acceleration  (Az)  along  the  z-axis  goes  through  a  similar  T/6 
sampling  and  averaging  process  and  is  feedback  at  a  T/2  sampling  interval 
through  digital  filters  (ITF)t/2,  (FTF)T/^»  and  (ATF)T/2.  This  same  pitch 
channel  loop  structure  is  used  for  the  yaw  rate  (r)  and  y-axis  linear  accel¬ 
eration  (Ay)  in  the  yaw  autopilot  loops.' 

C.  ROLL  AXIS  AUTOPILOT 

In  the  roll  axis,  roll  rate  (p)  about  the  missile  x-axis  is  also 
sampled  at  a  T/6  interval,  averaged  over  three  sanples  (filter  (JTF)t/6), 
and  feedback  at  a  T/2  interval  through  digital  filter  (NTTF)’1’/2.  The  roll 
angle  is  feedback  at  a  T/2  sampling  interval  through  digital  filter 
(LTF)t/2.  The  roll  autopilot  provides  effective  aileron  ccrtnands  (  6  a)  to 
the  actuators. 

D.  ANTI-ALIASING  LOW-PASS  FILTERING 

The  analog  signals  frcm  the  body  mounted  sensors  are  first  applied  to 
anti-aliasing  lew-pass  analog  filters  to  prevent  folding  of  the  high- 
frequency  sensor  noise  back  into  the  primary  frequency  band  of  the  analog- 
to-digital  (A/D)  samplers.  The  sensor  dynamics  and  anti-aliasing  filters 
are  treated  as  unity  transfer  functions  in  the  analysis  and  do  not  appear 
in  Figure  11. 


E.  SENSOR  THREE-SAMPLE  SIGNAL  AVERAGING 


The  output  of  each  anti-aliasing  filter  is  sanpled  at  a  T/6  interval. 
These  digital  signals  are  then  averaged  over  three  sanples  and  feedback  at 
a  T/2  interval  through  the  loop  compensation  to  the  servo  actuators.  This 
three  sanple  filtering  action  is  modeled  with  digital  filter  (JTF)t/6. 

In  the  discrete  time  domain,  this  filtering  action  is  given  by 

C(kT/6)  *  (l/3)R(kT/6)  +  (l/3)R[ (k-l)T/6]  +  (l/3)RC (k-2)T/6]  (72) 

In  this  difference  equation,  C(kT/6)  is  the  present  discrete  output  from 
the  filter,  and  R(kT/6)  is  the  present  sanpled  input  to  the  filter. 
R[(k-l)T/6]  and  RE  (k-2)T/6]  represent  the  sanpled  values  of  the  input 
obtained  during  the  previous  two  T/6  sampling  times.  This  difference 
equation  indicates  that  the  present  discrete  output  is  equal  to  the  numeri¬ 
cal  average  of  the  present  input  and  the  previous  two  sanpled  input  values . 

The  discrete  time-domain  difference  equation  can  be  modeled  directly  in 
the  z-plane.  The  exponent  on  the  z-^  delay  operator  is  interpreted  as  a 
time  domain  ordering  variable  for  each  discrete  signal  value.  The  coef¬ 
ficient  for  each  zT*  term  then  equals  the  value  of  the  discrete  signal  at 
the  kth  time  delay.  This  is  the  classical  sanpled  data  approach.  To  apply 
this  approach,  we  sinply  associate  each  discrete  time-domain  signed  value 
with  an  appropriate  z“*  delay  operator  (i.e.,  sinple  substitution). 

Equation  72  can  then  be  expressed  as 

C(z)  =  (l/3)R(z)  +  (l/Slz’Vz)  +  (1/3)z"2R(z)  (73) 

For  a  T/6  sampling  interval,  the  z- transform  variable  in  Equation  73  is 
defined  as  z  *  esT/6,  and  z T*  represents  a  one  sanple  period  delay  and 
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z“2  a  t>»  sanple  period  delay.  Hie  z”1r(z)  term  then  represents  the  value 
of  the  input  obtained  during  the  previous  sampling  time  and  the  z-2R(z) 
term  the  sampled  input  obtained  tvo  sanple  periods  earlier.  For  further 
insight  into  the  meaning  of  the  z“*  delay  operator,  see  Sections  II  and 
III. 

F.  DIGITAL  DATA  STORAGE  MODELING  VIA  ZERO-ORDER-HOLD 

Since  data  in  a  multi-rate  autopilot  is  updated  at  different  rates  and 
combined,  the  output  of  the  filters  in  the  slower  loops  must  be  stored  in 
the  microprocessor  software  between  filter  update  times.  For  example,  in 
Figure  11,  the  Azc  canmand  input  is  sampled  and  filtered  at  a  T  update 
interval  (i.e.,  every  T  seconds).  This  T  digital  signal  from  filter 
(KTF)t  is  combined  with  the  T/2  digital  signed,  from  filter  (ITF)t/2  in  the 
Az  acceleration  feedback  loop.  Since  the  T/2  signed,  arrives  at  the 
sunming  junction  twice  as  often  as  the  T  sampled  and  filtered  signal,  the 
latest  value  of  the  T  signal  from  filter  (KTF)t  must  be  stored  in  the 
software  between  the  T  update  intervals.  The  alternative  to  this  arrange¬ 
ment  is  to  operate  the  Azc  ocmnand  input  in  an  open- loop  mode  every  other 
T/2  update  of  filter  (ITF)T/2. 

In  the  multi-rate  autopilot,  the  digital  storage  of  data  between  updates 
is  modeled  in  the  s-plane  with  a  Zero-Order-Hold  (ZOH)  transfer  function, 


It  is  important  to  recognize  that  is  not  a  physical  ZOH  hardware  device 
(i.e.,  a  digital-to- analog  convertor).  It  simply  models  the  storage  of 
variable  data  within  the  microprocessor  software  between  the  algorithm 


update  times.  However,  the  M2  data  holds  in  Figure  11  represent  the  phys¬ 
ical  ZCH  hardware  units  that  implement  the  digital- to- analog  conversion 
of  the  signals  catmanding  the  actuators.  The  form  of  the  transfer 
function  is  the  same  as  the  functional  element. 


G.  COMMAND  AUGMENTATION  PREFILTER 

(KTF)t  in  Figure  12  is  a  command  augmentation  prefilter  that  is  used 
to  improve  the  steady-state  gain  and  dynamic  rerponse  of  the  missile  to 
acceleration  ocmnands  (Azc) .  Since  this  filter  is  outside  the  inner  pitch 
loop  and  outer  acceleration  loop,  it  has  no  effect  on  basic  missile  sta¬ 
bility. 

The  stability  margins  can  be  determined  by  opening  up  the  inner  loop  at 
point  E^/2  in  Figure  12  and  calculating  the  open-loop  z-plane  transfer 
function.  A  discrete  bode  plot  of  this  transfer  function  will  reveal  the 
gain  and  phase  margins. 


SECTION  V 


MULTI-RATE  MODELING  METHODOLOGY 


A.  INTRODUCTION 

Hie  multi-rate  modeling  methodology  presented  in  this  section  deals 
with  the  pitch  channel  autopilot.  Since  it  is  assumed  that  the  pitch  and 
yaw  autopilot  loops  are  identical  ( symmetrical  airframe  dynamics),  the 
results  in  this  section  also  apply  to  the  yaw  channel.  Hie  multi-rate 
equations  (transfer  functions)  representing  the  pitch  autopilot  are  for¬ 
mulated  directly  from  the  block  diagram  -in  Fi  gure  12 .  These  multi-rate 
equations  are  obtained  using  multi-rate  block  diagram  algebra. 

The  z-plane  closed-loop  and  open- loop  transfer  functions  for  the  pitch 
autopilot  are  developed  in  subsections  C  and  D,  respectively.  Discrete 
frequency  and  time-domain  performance  of  the  autopilot  can  be  assessed  from 
these  transfer  functions .  The  discrete  closed-loop  difference  equation 
follows  directly  from  the  closed- loop  z-plane  transfer  function.  This  dif¬ 
ference  equation  quantifies  the  time  domain  performance  of  the  autopilot 
(e.g.,  effective  time  constant,  rise  time,  peak  time,  settling  time,  and 
percent  of  overshoot).  To  complete  the  analysis,  closed- loop  bandwidth, 
steady  state  gain,  and  stability  margins  can  be  determined  by  substituting 
z  =  eT^T  in  the  z-plane  transfer  functions  and  calculating  the  discrete 
Bode  plots. 

B.  z-PLANE  PITCH  AUTOPILOT  BLOCK  DIAGRAM 

In  the  block  diagram  of  the  pitch  autopilot.  Figure  12,  the  third- 
order  actuator  model  ACT(s)  and  second-order  airframe  transfer  functions 


q(s)/DE(s)  and  Az(s)/q(s)  are  expressed  in  the  s-plane.  A  representative 
actuator  model  is  given  by 


ACT(s) 


_ Ka(wf)p _ 

(s  +  p)(s2  +  2  fas  +  w2) 

a  a  cl 


(76) 


where  ja  is  the  danping  ratio  and  wa  the  undanped  natural  frequency  of  the 
actuator.  For  missile  applications,  (s  +  p)  is  generally  a  large  pole 
(e.g. ,  p  >  100).  The  airframe  transfer  functions  would  appear  as 


q(s)  _  Kg  (s  +  a) 
^eTs) 


(77) 


Az(s)  _  Kz  (s  t  b)(s  -f  c) 
q(s)  Kq  (s+  a) 


In  Equation  77,  ^and  w  represent  the  short  period  danping  ratio  and 
undanped  natural  frequency  of  the  airframe,  respectively.  Az(s)  is  the 
acceleration  along  the  z-axis ,  q(s)  the  pitch  angular  rate,  and  DE(s)  the 
effective  elevator  deflection.  The  acceleration  transfer  function 
includes  a  ncnminimun  phase  zero  (i.e.,  a  zero  in  the  right-half  s-plane) . 

The  parameters  for  the  three  digital  ocnpensation  filters,  (ETF)T/2r 
(ETF)T/2f  and  (rrF)T/2f 

are  typically  scheduled  with  respect  to  the  dynamic 
state  of  the  missile  and/or  the  aerodynamic  environment  (e.g.,  Mach  number, 
altitude,  missile  velocity,  dynamic  pressure) .  Representative  forms  for 
these  z-plane  filters  can  be  expressed  as 

(DTF)T/2  .  +  aJ  (79) 

z2  +  a3 

r^T/2  _  ^  +  *2  tor\\ 


(81) 


(ITF)T/2  = 


Cjz2  +  c2 
z2  ^3 


In  Equations  79-81,  Z2  —  e®T/2  and  the  a's,  b's,  and  c's  are  numerical 
coefficients . 


The  digital  notch  filter  (ATF)'i'/2  is  used  to  suppress  the  first  body 
bending  mode  response  of  the  missile.  (ATF)T/2  generally  has  a  fixed  struc¬ 
ture  throughout  the  flight  envelope.  This  second-order  z-plane  filter  might 
appear  as 

(ATf)T/2  =  dkZ|  1  (82) 

z2  +  d4z2  +  d5 


Where  Z2  =  esT/2  and  the  remaining  variables  are  numerical  coefficients. 

The  three-sample  averaging  filter  (JTF)T/6  is  updated  every  T/6 
seconds.  The  output  of  this  filter  is  sampled  every  t/2  seconds  and  feed¬ 
back  through  the  compensation  in  each  loop.  The  fixed  structure  of  this  z- 
plane  filter  is  given  by 

2  . 

(jTF)T/6  =  - *6 - i  ,  =  eST/6  {83) 

3Z6 


Subsection  IV-E  describes  the  function  of  this  digital  filter. 


(KTF)t  is  a  command  augmentation  prefilter  that  is  used  to  improve  the 
steady-state  gain  and  dynamic  response  of  the  missile  to  acceleration  com¬ 
mands  (see  subsection  IV-6) .  The  general  form  of  this  filter  is 


tv _  elZ  *  ^  ,  «ST 


tQA  \ 


The  numerical  coefficients  for  the  (KTF)t  filter  are  normally  scheduled 
with  respect  to  the  dynamic  state  and/or  aerodynamic  environment  of  the 
missile. 


and  M2  in  Figure  12  are  Zero-Order-Hold  (ZOH)  s-plane  transfer  func¬ 
tions  expressed  as 


Mi  = 


(85) 


M^  models  the  process  in  which  the  latest  digital  output  from  filter 
(RTF)'1’  is  stored  in  the  microprocessor  software  between  the  T  update  times. 
The  T  sampled  signal  from  (KTF)t  is  combined  with  the  T/2  signal  from  the 
outer  acceleration  loop  to  farm  the  T/2  input  for  filter  (F*TF)T/2.  As 
pointed  out  in  Subsection  IV-F,  Mp  is  not  a  ZOH  hardware  element  in  the 
autopilot.  It  simply  models  a  storage  process  within  the  software.  M2 ,  on 
the  other  hand,  is  a  conventional  ZOH  hardware  device  that  converts  the 
digital  signals  from  the  microprocessor  to  analog  signals  for  the  actuators. 

C.  z— PLANE  CLOSED-LOOP  TRANSFER  FUNCTION 

This  subsection  sunmarizes  the  steps  necessary  to  form  the  multi-rate 
closed- loop  transfer  function,  Az/Azc,  for  the  pitch  autopilot  in  Figure  12. 
This  derivation  is  based  on  multi- rate  block  diagram  algebra.  For  the 
reader  who  is  unfamiliar  with  this  technique,  Appendix  A  presents  a  more 
detailed  derivation  of  the  following  multi-rate  equations.  The  following 
general  rules  of  thurib  are  utilized  in  this  subsection  and  Appendix  A: 

•  Start  at  the  inner  loop  and  work  outward. 


•  Designate  the  input  to  the  inner  loop  data  hold  as  a  discrete 


variable  (E^/2  in  Figure  12) 

•  Designate  the  input  to  the  inner  loop  sunning  junction  as  a  discrete 
variable  (E^/2  in  Figure  12) 

T 

•  Designate  the  sampled  ccninand  input  as  a  discrete  variable  (Azc  in 
Figure  12) 

Invoking  these  general  rules  will  reduce  the  computations  required  to  form 
the  closed- loop  multi- rate  transfer  function. 

We  start  the  derivation  by  writing  the  output  equation  for  the  inner 
loop  continuous  variable  g(s)  in  Figure  12. 

q(s)  =^fl)AOT(s)  M^s)  E?/2(z2)  (86) 

where 

Z2  *  e®T/2  (87) 

The  T/2  superscript  on  the  E^/2  discrete  variable  designates  the  sampling 
interval  in  seconds.  The  update  rate  far  each  filter  and  sample  rate  for 
each  signal  within  the  autopilot  block  diagram  will  be  designated  using  this 
multi-rate  superscript  notation  (see  Section  III).  In  addition,  for  simplic¬ 
ity,  the  (s)  and  general  (z)  notation  for  each  transfer  function  will  be 
dropped  far  the  remaining  derivation.  For  example,  the  filter  ATF( Z2 )  is 
designated  (ATF)'r/2,  vfriere  T/2  is  the  sanpling  or  update  interval  as 
reflected  in  the  z- transform  variable  Z2  =  esT/2. 

Utilizing  multi- rate  block  diagram  algebra,  the  loop  equation  for  the 
discrete  variable  E?/2  can  be  expressed  as 


E?/2  =  (ATF)t/2e£/2  +  (ATF)T/2(DTF)T/2  [(JTF)T/6(-^g  ACT  M2)T/6]T,/2Efr/2 


=  [1  -  (ATF)T/2(DTF)T/2[(JTF)'I’/6(-2g  ACT  M2)T/6]T,/2]  1  (ATF)T/2eT/2 


=  GT/2^/2 


(88) 


with 

qT/2  =  [1  -  (ATF)T/2(DTF)T/2C(JTF)T/6(-^  ACT  M2)T/6]T//2]  1  (ATF)t/2  (89) 


Substituting  Equation  88  into  Equation  86  gives  the  closed- loop  transfer 
function  far  the  inner  pitch  rate  loop  (Equation  90) . 

q  =  (^g  ACT  M2)gT/2eT/2  (90) 

Then 

q  =  G*Ej/2  (91) 

with 

G*  =  (|g  ACT  M2)Gj/2  (92) 


In  Equation  90,  the  first  bracketed  tern  on  the  right-hand  side  is  a  s-plane 
transfer  function  and  G^/2  is  a  z- plane  transfer  function.  Expressing 
Equation  90  with  more  explicit  notation  gives 

a&b =  OT(S)  "2(s)  <5/2(Z2)  ,93) 

1  ^ 

We  next  write  the  outer  loop  equation  for  the  output  continuous  variable 
A z,  with  the  inner-pitch  rate  loop  closed  (Figure  13) .  This  outer  loop 
equation  is  given  by 


hz  =  (££>  G*eT/2 


(94) 


ms 


As  in  the  inner  loop,  the  discrete  variable  E^/2  can  be  expressed  as 


ET/2  _  (PTF)T/2Mj/2(KTF)TAzc  -  (PrF)T/2(lTF)T/2[(JTF)T/6(^  G*)t/6]T/2e]'/2 
=  [1  +  (PTF)t/2(ITF)t/2[(JTF)t/6(^  G*)T/«3T/2]  1  (Fir)T/2^/2(KTF)TAZc 

«  GT/^T&ic  (95) 

with 

gT/2  =  [1  +  (FTF)t/2(ITF)t/2[(JTF)t/6(-^  G*)t/6]T//2]  1  (PTF)’1’/2^’/2  (96) 

G£  =  (KTF)T  (97) 


The  overall  closed- loop  transfer  function  is  obtained  by  substituting 
Equation  95  into  Equation  94. 

M  g*)gT/2gJ  (98) 

Reverting  back  to  a  more  explicit  notation,  Equation  98  becomes 

7^r7T  =  qrar'  ra?sT,Cr<S>  M2<s)  Ga/2(z2)GB/2(z2>  GC(z)  (99) 

AZcl S ) 

and 

z  =  esT  ,  Z2  =  esT/2  (100) 

Taking  the  z- trains  form  of  both  sides  of  Equation  99  with  respect  to  a  T 
sampling  interval  produces 

§E  *CT  ”2  °C  (101) 
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Notice  in  Equation  101  that  the  bracketed  factor  is  a  high- to- lew  rate  trans¬ 
form  conversion  and  G^/2  and  g^/2  cannot  be  separated  frem  the  T  resampling 
operation  (see  Section  III).  That  is, 

f  k  **  "2  ^/2gI/2]t  »*  f  b  m  »2>t  gI/2gI/2  «°2> 

Employing  the  phantom  sampler  principle  in  Subsection  III-E,  Equation  101 
can  be  expressed  as 

-  C  <f  ^  act  Mj )  T/W/M^T  GT  ( 103 ) 

Azc 


where 

GJ/2  ~  ^  "  (atf)T/2(dtf)T//2  C(JTF)t/6  (-^  ACT  M2)T^6]T^2^  1  (ATF)t/2  (104) 

G^/2  =  [1  +  (FTF)T/2(ITF)T/2G^/2  [(JTF)T/6  (•—  ACT  M2)T/S]'r/23  1  (PrP)T/2^/2 

(105) 

Gj  =  (KTF)t  (106) 

Equation  103  is  the  z-plane  closed- loop  transfer  function  for  the  pitch  auto¬ 
pilot  with  z  defined  for  a  l/T  sample  rate  (i.e.,  z  =  esT) .  A  nunerical 
exanple  that  illustrates  the  evaluation  of  Equation  103  in  great  detail 
appears  in  Appendix  B.  Also,  an  automated  technique  for  evaluating  Equation 
103  is  presented  in  Appendix  C.  This  automated  evaluation  utilizes  the 
macro  capability  of  the  linear  analysis  program  TOTAL  (Reference  13). 

A  sanpled  step  input 

=  <i>T  -  TTT  (107) 
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can  be  applied  to  Equation  103  to  determine  the  discrete  time  domain 
figures  of  merit  for  the  pitch  autopilot.  Discrete  time  response  tech¬ 
niques  are  discussed  in  subsection  II-F. 


The  sampled  spectrun  closed- loop  frequency  response  is  obtained  by  eval¬ 
uating  Equation  103  at 


z  =  =  1  ^£wT  (108) 

This  evaluation  is  expressed  as 

z  =  eST 

z  =  i2£wr  (109) 


A  +  jB  -C(i£  ACT  t*)T/2dT/J(£/S]T  gt 


-  Mag  y  <6 


where 


Mag  =  (A2  +  B 2)^  ,  jb  =  Tan_1(B/A)  (110) 

A  _ 

The  z  *  e31  notation  in  Equation  109  is  for  the  purpose  of  calling  out  the 
definition  of  z  being  used  in  the  evaluation.  The  complex  z-plane  evaluation 
in  Equation  109  at  z  *  eJwT  is  analogous  to  the  classical  evaluation  of  s  =  j 
in  the  s-plane.  That  is,  at  any  given  sinusoidal  input  with  frequency  w  (see 
Table  1) , 


(sin  wT)t  = 


_ z  sin  v/T _ 

z2  -  (2  cos  wT)z  +  1 


z  =  e®^ 


(111) 


the  output  is  also  sinusoidal,  with  frequency  w,  differing  from  the  input 
by  only  a  magnitude  factor  (Mag)  and  phase  angle  (  0) .  This  represents 

the  traditional  concept  of  the  '  sampled  spectrum'  which  is  limited  to 
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determining  the  single  sinusoid  that  fits  the  system  output  at  the  sample 
points.  For  a  more  ccmpr ehens i ve  concept  of  multi-rate  frequency  response 
see  References  1,  11,  12,  14,  and  15.  The  numerical  frequency  evaluation 
given  by  Equation  109  is  implemented  as  Options  50-59  and  150-159  in  the 
linear  analysis  ccmputer  program  TOTAL  (Reference  13). 

D.  Z-PIANE  OPEN-LOOP  TRANSFER  FUNCTION 

The  open- loop  transfer  function  for  the  pitch  autopilot  is  calculated 
in  this  subsection  with  the  loop  break  point  at  E^/2  in  Figure  12 
(sometimes  referred  to  as  the  servo  break  open- loop  transfer  function). 

At  this  break  point,  two  loop  paths  define  the  open- loop  z-plane  transfer 
function.  These  two  loops  are  given  by 

(OLTF)t/2  =  (ATF)T/2(ETF)T/2(lTF)T/2C(JTF)'r/6(^  ACT  M2)T/6]T/2 

-  (ATT)T/2(DITOT/2[{.riF)T/6(^  ACT  m^T/G-jT /2  (n2) 


The  negative  sign  far  the  second  term  (loop)  in  Equation  112  is  a  result 
of  the  positive  sign  on  the  inner  loop  sunning  junction. 

The  sanpled  spectrum  open- loop  frequency  response  is  obtained  by  eval¬ 
uating  Equation  112  at  z  =  e3wT  =  1<£>T.  This  evaluation  is  expressed  as 


A  +  jB  =  Gj/2  -  G^/2 


=  Mag^0 


z  i  e^/2 
z  =  I  wT 


(113) 


or 
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Z  4  eST/2 


A  +  jB  =  G^/2 

z  =  l2£wT  (114) 

=  Mag  2^.0 

where 

Mag  =  (A2  +  B2)3^  ,  =  Tan-1(B/A)  (115) 

and 

G^/2  =  (ATF)T/2(FrP)'T’/2(rrF)T/2  [(JTF)T/6(-g|  act  M2)t/6]T/2  (116) 

gT/2  =  (ATF)t/2(DTF)t/2C(JTF)t/6(^  ACT  M2)T/6]T/2  (117) 

dT/2  =  GT/2  -  g|/2  (118) 

When  the  magnitude  (Mag)  and  phase  (0)  frcm  Equation  113  or  114  are 
plotted  with  respect  to  frequency  (w) ,  the  discrete  Bode  plots  are  obtained. 
These  plots  can  be  used  to  assess  the  gain  and  phase  margins  for  the  pitch 
autopilot.  In  the  TOTAL  ccrrputer  program,  Cptions  50-59  and  150-159  can 
be  used  to  evaluate  Equations  113  and  114. 


SECTION  VI 


COMPUTATIONAL  ALGORITHMS  FOR  MULTI-RATE  SYSTEMS 

A.  INTRODUCTION 

The  basic  problem  facing  the  control  engineer  is  to  obtain  valid 
discrete  models  of  multi-rate  hybrid  systems  (i.e. ,  systems  containing  both 
analog  and  multi-rate  discrete  elements).  These  multi- rate  models  must  be 
in  a  convenient  form  that  can  be  readily  analyzed  using  the  analysis  and 
synthesis  tools  available  today.  The  preceding  sections  of  the  report  pre¬ 
sented  a  methodology  far  feruling  these  multi- rate  models.  The  following 
subsections  discuss  computational  algorithms  that  are  well  suited  for 
numerically  evaluating  these  models.  These  algorithms  are  also  utilized  in 
the  evaluation  technique  presented  in  Appendix  C. 

B.  DISCRET  COMPUTER  PROGRAM 

Once  a  valid  multi-rate  model  equation  is  obtained,  the  next  fundamen¬ 
tal  step  is  to  evaluate  the  model  and  place  it  in  a  convenient  transfer 
function  farm.  This  includes  the  discretization  of  the  continuous  elements 
in  the  system  into  a  valid  discrete  domain.  The  DISCRET  computer  program 
presented  in  Reference  11  provides  this  discretization.  DISCRET1  takes  a 
continuous  element  expressed  as  an  s-plane  transfer  function  and  transforms 
it  into  the  z- ,  w- ,  or  w' -plane.  The  program  can  calculate  the  standard, 
delayed,  or  advanced  discrete  transform.  Data  holds  including  the  zero 
order,  first  order,  second  order,  and  slewer  can  be  inserted  into  the  trans¬ 
formation. 


The  discrete  models  of  the  continuous  elements  in  the  multi-rate  model 


equation  can  be  readily  combined  with  the  inherent  discrete  elements  (e.g. , 
digital  filters)  to  provide  a  unified,  concise  description  of  the  system. 

An  inherently  discrete  element  may  first  be  modeled  with  a  difference 
equation  and  then  directly  converted  to  the  z-plane  by  substituting  the 
z“k  delay  operator  for  each  discrete  term  in  the  difference  equation. 
Therefore,  all  the  elements  in  the  multi-rate  model  equation  can  be 
described  by  transfer  functions  in  the  z-,  w-,  or  w' -plane.  Consequently, 
similar  frequency  and  time  domain  techniques  used  for  continuous  s-plane 
systems  can  be  applied  by  the  control  system  engineer  to  a  wide  variety  of 
practical  multi-rate  systems. 

C.  TXOCNV  COMPUTER  PROGRAM 

In  the  nulti-rate  model  equations,  a  troublesome  situation  arises  when 
the  digital  signals  are  resampled  at  a  lower  rate.  This  sampling  con¬ 
figuration  requires  the  calculation  of  a  low-rate  z-plane  transform  ffan  a 
given  high-rate  z-plane  transfer  function.  The  TXCONV  computer  program  pre¬ 
sented  in  Reference  11  calculates  this  ccrtplex  transformation  and  makes  it  a 
routine  operation  in  the  analysis  and  design  process.  Specifically,  TXOONV 
takes  a  high-rate  discrete  transfer  function  expressed  in  the  z-,  w-,  cr 
w' -plane  and  transforms  it  into  a  desired  lew-rate  discrete  transfer  func¬ 
tion  in  the  z- ,  w-,  or  w'  -plane. 

The  fundamental  definition  of  the  z- transform  and  the  discrete  inver¬ 
sion  integral  (see  Section  II)  form  the  basis  for  the  algorithms  in  the 
TXCONV  computer  program.  The  computer  mechanization  of  the  high- to- lew  rate 
transform  conversion  relies  on  the  practiced,  calculation  of  the  residues 
of  a  ocmplex  integral.  The  residues  for  this  integral  are  calculated  in  an 
unconventional  manner  using  a  limiting  process  via  L'Hopital's  rule.  This 


method  simplifies  the  mechanization  scheme  and  leads  to  a  closed- form 
solution. 

D.  TOTAL  COMPUTER  PROGRAM 

Both  the  DISCRET  and  TXOONV  computer  programs  (Figure  14)  represent 
essential  tools  far  dealing  with  multi-rate,  hybrid  control  systems.  They 
provide  two  of  the  primary  ccrputational  algorithms  necessary  to  evaluate  the 
multi-rate  model  equation.  Hcwever,  they  are  not  the  only  tools  or  tech¬ 
niques  required  by  the  analyst.  Other  ccrputational  routines  that  are 
required  include  those  that  can  manipulate  and  combine  the  polynomials  and 
transfer  functions  in  the  multi-rate  model  equation.  In  addition,  analysis 
algorithms  are  needed  to  calculate  the  root  locus,  frequency  response,  and 
time  response  in  the  z-,  w-,  or  w1 -plane.  All  these  computational 
algorithms  are  available  in  the  TOTAL  computer  program  (Reference  13). 

DISCRET  and  TXOONV  have  also  been  integrated  into  TOTAL. 

Over  100  options  are  new  available  to  the  user  during  the  interactive 
execution  of  the  TOTAL  program.  This  interactive  feature  allows  a  close 
coupling  between  the  analyst  and  the  computing  machine  (digital  computer) , 
and  a  real-time  dialog  between  the  two  can  be  effectively  carried  out.  This 
results  in  a  more  effective  and  efficient  usage  of  the  computing  machine  and 
improves  the  accuracy  and  speed  of  the  analysis  and  synthesis  process. 
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Figure  14.  General  Structure  of  DISCRET  and  TXCONV 
Computer  Programs 
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APPENDIX  A 


EXPANDED  MULTI-RATE  MODELING  METHODOLOGY 


INTRODUCTION 

This  appendix  expands  the  derivation  presented  in  Section  V  that  leads 
to  Equation  103.  The  steps  in  the  following  derivation  are  listed  in  great 
detail.  Before  proceeding,  it  will  be  helpful  to  review  the  lcw-to-high  and 
high- to- lew  rate  sampling  identities.  The  general  lcw-to-high  rate  sampling 
scheme  is  represented  as  (N  is  a  positive  integer) 

C(^(z)]T/N  =  Gt(z)  '  ,  z  =  e57  (A-l) 

However,  the  reverse  of  Equation  A-l  is  not  true.  That  is, 

[Gt/n(zn)]T  +  gT/N(Zn)  ,  zn  =  eST/N  (A-2) 

Equation  A- 2  represents  high- to- low  rate  sampling  in  which  a  high-rate 
sargpled  function  Gt/n(zjj)  is  resanpled  at  a  lower  rate  (l/T  Hz).  These  two 
identities  can  be  further  expanded  into  the  following  identities: 

CG1(s)GJ(z)]T/N  =  gT/N(Zn)gJ(z)  (A-3) 

[gT/N(Zn)GT(z)]  t  =  CgT/N(zn)]  T  GT(z)  (A— 4) 

DETAILED  DERIVATION 

Figure  12  from  Section  V  is  repeated  here  as  Figure  A-l.  The  update 
intervals  (i.e.,  T,  T/2,  and  T/6)  are  shown  with  each  digital  filter.  The 
intent  of  the  following  derivation  is  to  obtain  the  closed- loop  z-plane 
transfer  function  given  by 


(A-5) 


if  -  C(£§  KT  M2)T/2GT/2GT/2f  qT 
Azc 


where 

GT/2  =  [1  -  (ATF)T/2(DTF)T/2[(JTF)T/6(|L  ACT  )T/6]T/2]  1  (ATF)T/2  (A-6) 

gT/2  =  [1  +  (FTF)T/2(ITF)T/2G'T/2[(JTF)T/6 (^|  ACT  M 2)T/e]T/2]  1  (FTF)T/2MT/2 

(A-7) 

gT  =  (KTF)t  '  (A-8) 

We  start  the  derivation  at  the  input  to  the  M2  data  hold  in 
Figure  A-l.  Hie  analog  pitch  rate  can  be  expressed  as 

q  -  ijL  ACT  M2  E^/2  (A— 9) 

Equation  A-9  indicates  that  the  discrete  signal  ET/2  is  multiplied  by  the 
analog  transfer  functions  M2*  ACT,  and  q/DE  in  order.  M2  is  a 
Zero-Order-Hold  (ZOH)  that  converts  the  discrete  signal  ET/2  into  an  analog 
signal  for  the  actuator,  ACT.  The  airframe  transfer  function  q/DE  relates 
pitch  rate  to  elevator  deflection. 

The  discrete  signal  ET/2  can  be  expressed  as  the  son  of  two  signal 
paths,  with  one  signal  path  being  a  complete  loop.  One  signal  path  is  asso¬ 
ciated  with  the  discrete  signal  E^/2  and  the  other  with  the  ET/2  signal 
itself.  For  the  first  signal  path,  the  E^1/2  signal  is  simply  multiplied  by 
the  (ATF)t/2  digital  filter. 


In  the  second  signal  path  which  forms  a  complete  loop,  the  ET/2 

signal 

is  first  multiplied  by  M2 ,  ACT,  and  q/DE. 

3g  ACT  M2e'1'/2  (A— 11) 

The  product  in  Equation  A-ll  is  next  sampled  every  T/6  seconds  by  sampler 
No.  1. 

ACT  M2Et/2]T/6  (A-12) 

Utilizing  Equation  A-3 ,  Equation  A-12  reduces  to 

[■2g  ACT  M2  Efr/2]T^6  =  ACT  M23t/6  Et/2  (A-13) 

The  product  in  Equation  A-13  is  then  multiplied  by  digital  filter 
(JTF)t/6  and  resampled  at  a  T/2  interval  with  sampler  No.  2.  These 
operations  are  expressed  as 

[(JTF)t/6(|L  act  M2)t/6Et/2]T//2  (A-14) 

The  lower  rate  sampled  signal  E?/2  can  be  pulled  out  of  the  bracket  in 
Equation  A-14  (see  Equation  A-4) .  However,  since  we  are  resampling  the 
remaining  high-rate  factors  at  a  lower  rate  (i.e.,  at  2/T  Hz),  Equation 
A-14  becomes 

[(JTF)T/6(-2g  ACT  M2)t/6Et/2]T/2  =  [(jTF)T/6(-2g  ACT  M2)t/6]T/2  E?/2  (A-15) 

When  Equation  A-15  is  evaluated,  the  high- to- lew  rate  transform  will  be 
numerically  calculated.  For  now,  we  simply  carry  along  this  bracketed  fac¬ 
tor.  To  complete  the  loop,  the  product  in  Equation  A-15  is  multiplied  by 


digital  filters  (OTF)t/2  and  (ATF)T/2. 

(ATF)T/2(DTF)T/2C(JTF)T/6(-3g  ACT  M2)T/6]T/2  E?/2  (A-16) 

Sunming  Equation  A-10  and  Equation  A-16,  the  E^/2  sampled  signal  is  given  by 

E?/2  =  (ATF)t/2  E^/2  +  (ATF)T/2)(DTF)'r/2C(JTF)'r/6(-^  ACT  M2)T/6D'r^2  E?/2 

(A-17) 

To  solve  far  E^/2,  we  first  collect  terms  on  the  left  hand  side  of  the 
equation. 

[i  -  (\tf)t/2(dtf)t/2[(jtf)t/6(-^  act  m2)t/6]t/2]et/2  =  (atf)'i’/2eJ/2  (A— 18) 


Multiplying  both  sides  of  Equation  A-18  by  the  inverse  produces 

EfT/2  =  [1  -  (ATF)T/2(DTF)T/2[(JTF)T/6(^  ACT  M2)T/6:T/2]_1  (ATF)T/2ri[/2 

(A— 19) 


Let 

gT/2  =  [1  -  (ATF)t/2(DTF)t/2[(JTF)t/6(-^  ACT  M2)T/6]T//2]  1  (ATF)t/2  (A-20) 


Then 

EfT/2  =  GT/2ET/2 


Substituting  Equation  A-21  into  Equation  A-9  gives 

q  =  ACT  Ms  Gj/2^/2 


(A-21) 


(A-22) 


If  we  next  let 


G*  -  h  «2  <S/2 


(A-23) 


l 


Equation  A-22  then  becomes 

q  =  G*E^/2  (A-24) 

As  pointed  out  in  Section  V,  G*  is  a  hybrid  function,  that  is,  it  contains 
both  analog  and  discrete  factors.  The  ACT  M2  product  results  in  a  s- 
plane  transfer  function  and  G^/2  is  a  z-plane  transfer  function  with 
z  =  esT/2.  Utilizing  Equation  A-24,  we  can  redraw  Figure  A-l  as  Figure  A-2. 

The  next  ..ajar  step  in  the  derivation  is  the  closing  of  the  outer  accel¬ 
eration  loop.  Fran  Figure  A-2,  the  analog  normal  acceleration  is  given  by 

Az  =  G*Ej/2  (A-2S) 

Similar  to  E^/2,  the  discrete  signal  E^/2  can  be  expressed  as  the  son  of  two 

signal  paths,  with  one  signal  path  being  a  complete  loop.  The  first  signal 

.  .  .  T 

path  is  associated  with  the  sampled  ccrmand  input  Azc  and  the  second  with 

E^*/2  itself.  The  first  signal  path  is  given  by 

(prF)T/2[Ml(KTF)TA|c]T/2  (A-26) 

In  Equation  A-26  the  Azc  sampled  signal  is  first  multiplied  by  (KTF)t  and 
.  This  product  is  resampled  at  a  t/2  interval  and  then  multiplied  by 
digital  filter  (ETF)t/2.  Since  this  resampling  is  at  a  higher  rate. 

Equation  A-26  reduces  to  (see  Equation  A-3) 

(FTF)t/2mT/2(KTF)tAZc  (A-27) 
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The  second  signal  path  associated  with  El^/2  is  obtained  by  tracing  com¬ 
pletely  around  the  outer  acceleration  loop.  The  E^/2  discrete  signal  is 
first  multiplied  by  the  G*  and  Az/q  transfer  functions. 

&  G*E?/2  (A-28) 

q  i 

This  product,  representing  the  signal  at  the  output  of  the  Az/q  block  in 
Figure  A-2 ,  is  next  sanpled  with  a  T/6  interval  (sampler  No.  3). 

[|E  G*eT/2]T/6  (A-29) 

In  Equation  A-29,  the  E^/2  signal  can  be' pulled  out  since  we  are  resampling 

at  a  higher  rate  (i.e.,  6/T  Hz). 

(Az  G*)T/6ET/2  (A- 30) 

q  i 

The  product  in  Equation  A-30  is  multiplied  by  digital  filter  (JTF)T/6  and 
the  results  is  resampled  every  T/2  seconds  (sampler  No.  4). 

[(JTF)T/6(^  G*)t/6eT/2]T/2  (A-31) 

Pulling  the  E^/2  signal  out.  Equation  A-31  becomes 


C(JTF)t/6(—  g*)t/6]T'/2ET/2  (A-32) 

q  i 

The  bracketed  factor  in  Equation  A-32  is  a  high- to- lew  rate  conversion 
which  must  be  numerically  calculated  (as  is  also  true  for  Equation  A-15). 

The  final  step  for  this  signal  path  is  the  multiplication  by  digital  filters 
(ITF)T/2  ^  (ftf)t/2. 
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(FIT-)T/2(iTf)T/2[{JTP)T/6(^g*)T/6]T/2eT/2 


(A-33) 


Canbining  the  two  signed  path  Equations  A-27  and  A-33  produces 


E^/2  =  (Fir)T/2^/2(KTF)'rA2c  -  (FrF)T/2(lTF)T/2C(jTF)T/6(^  g*)t/6]T/2e£/2 

(A-34) 

Solving  for  E^/2  gives 


E^/2  =  [1  +  (liTF)T/2(lTF)T/2[(JTF)'r/6(-^  G*)T/6]T/2]  1  (FTF)T/2MjY2(]^)TA^ 

(Ar-35) 


If  we  let 

G|/2  =  [1  +  (FTF)T/2(lTF)T/2[(JTF)T/6(^  G*)t/6]T/2]  1  (FTF)t/2^/2  (a-36) 


g£  =  (KTF)t 


(A-37) 


Then 

EiT/2  =  qT/Zg^ 


(A- 38) 


The  overall  closed  loop  transfer  function  is  obtained  by  substituting 
Equation  A-38  into  Equation  A- 2 5 . 

Az  =  M  G*gT/2gTa£c  (A-39) 

We  next  take  the  z- transform  of  both  sides  of  Equation  A-39  with  respect  to 
a  T  sanpling  interval  (z  =  esT) . 

Al=  G*gT/2gTA!c)T 


(A-40) 


The  T  transform  factors  can  be  pulled  out  and  Equation  A-40  becomes 


ijjj  -  <£  G*GT/2)TGT  (Ml) 

*Ie  q 

Substituting  for  G*  (Equation  A-23) ,  Equation  A-41  can  be  expressed  as 

4  -  W2>T<£  ««*) 

Azc 

The  final  step  in  the  derivation  is  to  introduce  a  phantom  T/2  sampler  to 
Equation  A-42.  This  mathematical  operation  is  depicted  in  Figure  A-3.  This 
step  is  valid  since  the  T  output  sampler  simply  rejects  all  the  unwanted 
samples  from  the  T/2  phantom  sampler.  Equation  A-42  then  becomes 

J  -  [<§  ACT  M2)T/2gT/20T/2]t  gT  (M) 


Equation  A-43  is  the  final  result,  the  closed-loop  z-plane  transfer  function 
far  Figure  A-l  (z  =  esT) .  The  intermediate  transfer  function  products  G^/2, 
(g/2,  and  gT  are  given  by  Equations  A-20,  A-36,  and  A-37,  respectively. 
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APPENDIX  B 


NUMERICAL  EVALUATION  OF  MULTI-RATE  MODEL  EQUATIONS 

INITODUCTION 

This  appendix  outlines  the  step-by-step  procedures  for  nurerically  eval 
uating  the  z-plane  closed- loop  transfer  function  expressed  in  Equations 
103-106,  Section  V.  This  multi-rate  equation  is  repeated  below. 

A|  -  C<|§  act  M2)T/2gT /2<q/2f  gt  (b-1) 

G^/2  =  [1  -  (ATF)t/2(DTF)t/2C(JTF)t/6(^  ACT  M2)T/6]T^2]  1(ATF)t/2  (B-2) 

Gg/2  =  Cl  +  (FrF)T/2(ITF)T/2Gj/2[(jTF)T/6(^  ACT  M2)t/6]T/2]’*1(FTF)t/2I^'/2 

(B-3) 

Gj  =  (KTF)t  (B—4) 

Equations  B-2  and  b-3  are  first  calculated  and  then  substituted  into 
Equation  B-l.  A  canplete  nunerical  exanple  is  presented  to  help  clarify 
each  step  in  the  evaluation  process  and  to  being  the  mathematical  details 
sharply  into  focus. 

MAJOR  STEPS  IN  EVALUATION 

The  following  transfer  functions  and  sanple  rates  are  used  in  the  sca¬ 
lar  illustrative  exanple: 

(ATF)  T/2  =  (ErTF)T/2  =  (FTF)T/2  =  T/4^z2  +  ^ 


T 


(lTF)T/2  =  (KTF)t  =1.0 

(B-6) 

(JII.,T/6  = 

(B-7) 
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q  _  1_  Az  _  1 

DE  8  '  DE  ~  ^5 


ACT  =  1.0  (B-8) 


Mj.  = 


(B-9) 


z2  *  eST/2  ,  zg  =  eST/6  (B-10) 

T  =  1/40  ,  T/2  =  1/80  ,  T/6  =  1/240  (B-ll) 

Hie  evaluation  of  Equations  B-l  through  B-4  can  be  separated  into  nine 
major  items.  Each  of  these  items  require  a  series  of  nvmerical  calculations . 
These  major  itens  are  listed  next. 


(1)  [(JTF)t/6(|L  ACT  M2)T/6]t/2 

(2)  (ATF)T/2(DrF)T/2C(JH?’)T/6(^  ACT  M^T/e^ 

(3)  Cl  -  (atf)t/2(dtf)t/2C(jtf)t/6(^  ACT  M2)T/6]T/2]-l 

(4)  G^/2  *  [1  -  (ATF)T/2(DrF)T/2[(JTF)T/6(-3L  ACT  M2)T/6;jt/2]  1(ATF)t/2 

(5)  [(jtf)t/6(£§  per  m2)t/6]t'/2 

(6)  (PrF)T/2(lTF)T/2Gj/2C(JTF)T/6(^  ACT 
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(7)  [1  +  (FTF)'r/2(n,F)T/2Gj/2C(JTF)T/6(-^  ACT  M2)T/6]T/'2] 

{B)  G^/2  «  [1  +  (FTF)T/2(ITF)T/2G^/2C(JTF,)T/6(^|  ACT  M2)T/6]T'/2]  1(FTF)t/2mT/2 


(9)  ^  -  [(M  jct  MjIT/^T /J^/2]T  gT 
AZc 

In  items  1,  5,  and  9,  a  high- rate  z-plane  transfer  function  must  be  con¬ 
verted  to  an  equivalent  lew-rate  z-plane  transfer  function.  For  example,  in 
item  1,  the  T/2  trains  form  is  obtained  from  a  T/6  transform.  That  is, 

C(JTF)T/6(S_  ACT  M2)T/6]T/2  =  [Gt/6(z6)]t/2  „  gT/2(z2)  (B-12) 

DE  1  2 


where 

=  esT/6 

and 

z2  =  esT/2 


(B-13) 


(B-14) 


A  number  of  techniques  are  available  for  calculating  the  standard 
Zg- transform  of  the  s-plane  function  q/DE  ACT  ^  [i.e.,  (q/DE  ACT  M2)T/6]. 
The  most  cannon  technique  is  to  apply  partial  fraction  expansion  to  the  s- 
plane  function  and  then  utilize  a  table  of  z- transform  (Table  1).  A  com¬ 
puter  program  (ADVANZ)  that  implements  this  technique  is  described  in 
Reference  11.  The  source  listing  for  this  program  can  be  found  in 
Reference  15  and  is  implemented  as  options  141-143  in  TOTAL  (Reference  13) . 
This  program  converts  a  general  s-plane  transfer  function  into  the  z-,  w-, 
or  w* -plane.  It  can  calculate  the  standard,  delayed,  or  advanced  discrete 
transform;  and  can  insert  a  zero  order,  first  order,  second  order,  or 


slewer  data  hold  into  the  transformation. 

The  theoretical  basis  (and  a  computer  program)  for  performing  the  high- 
to-lcw  rate  z-plane  transform  conversion  in  Equation  B-12  is  described  in 
Reference  11.  The  source  listing  for  this  computer  program  (TXCONV)  appears 
in  Reference  15  and  is  implemented  as  options  144-149  in  TOTAL.  This 
program  allows  a  high- rate  transform  in  the  z-,  w-,  or  w*  -plane  to  be  con¬ 
verted  to  a  lew- rate  discrete  transform.  The  user  simply  inserts  the  high- 
rate  transfer  function  G^/6  and  the  sampling  ratio  (T/2)/(T/6)  =  3  and  the 
program  calculates  the  lew-rate  transform  G^/2. 

NUMERICAL  CALCULATIONS 

In  this  subsection,  the  actual  numerical  calculations  for  items  1-9  are 
carried  out.  It  is  normally  necessary  to  utilize  linear  analysis  software 
aids  (e.g. ,  TOTAL)  to  perform  these  calculations.  However,  the  following 
calculations  will  be  carried  out  essentially  by  hand. 

ITEM  1  [(JTF)t/6(-^  ACT  M2)t/6]T/2 

Substituting  the  appropriate  transfer  functions  from  Equations  B-5  through 
B-9  produces 

(Sg  ACT  M2)t/6  =  (1  -  e  I)T/6  (b_15) 

Since  zg  =  eST/6, 

e-sT/2  =  (esT/6)-3  =  z-3  (b-16) 

Replacing  e-sT/2  in  Equation  B-15  with  z£3  and  pulling  out  the  resulting  z- 


plane  factor  gives 


(3-  ACT  M2)T/6  =  (1  -  z-63)(VT/6 


(B-17) 


In  Equation  B-17,  the  T/6  z- transform  of  the  s-plane  transfer  function 
l/s2  is  required.  Frcm  a  table  of  z- transforms  (Table  1), 


&  ACT  ^(T/6  _  (a  -  Zg3) 


(Z6  -l)(z|  +  Z6  +  1)  (T/6)ze 

4  <«6  -  n: 


(B-18) 


C(jtf)t/6(-3-  ACT  M2)t/6]T/2  = 


rzg  +  Zg  +  1  (T/6)(z|  +  zg  +  l)nT/2 

L - 5 -  - 5 - J 

3z6  _  15 


^(T/6) (zg  +  zg  +  l)2j 
3zg(zg  -  1) 


(B-19) 


Hie  T/2  superscript  notation  in  Equation  B-19  tells  us  to  take  the 
Z2  =  e87/2  transform  of  the  function  in  the  brackets.  Since  this 
function  is  a  high- rate  z-plane  transfer  function  (zg  =  esT/^),  a  high-to- 
low  transform  conversion  is  needed.  Fran  Section  III,  Equation  71,  this 
conversion  can  be  obtained  by  combining  the  residues  of  the  following 
expression. 


<%/*(  Zk)  =  [<^/P(zp)]T/k  =  ^residues 


Gi/p(zp)  Zk 


Zk  -  ZN 


(B-20) 


Zp  ss  poles  of 
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N  =  ratio  of  lcw-to-high  rate  sanpling  intervals  (p/k) 


G^/^Czj^)  =  lew- rate  z-plane  function  (z^  =  esT/^) 


<jj/P(«p)  =  high- rate  z-plane  function  (zp  =  esT/P) 


For  Equation  B-19, 


z^  =  z2  =  esT/2  ,  Zp  =  z£  =  esT/6  ,  N  =  3 


G^/p(Zp)  (JTF)T//6(jL  ACT  M2)T/6  (T/6)(z|  +  Zg  +  l)2  (B 

*p  3z|(z6  “ 


Substituting  Equation  B-21  into  Equation  B-20  gives 


[(jTF)T/6(^  act  M2)T/6]T/2  =Y~reg  (T/6)  (|L+  ^  + 
UC.  *  * -  -  1  ) 


1)  z2 


3z|(ze  -  1)  z2  -  z3 


ze  =  ° 


The  residue  calculations  can  be  sinplified  by  recognizing  that 

z|  =  z|z2  *  ZjZ^.  The  z2  in  the  numerator  of  the  z2/(z2  -  z3)  factor  can 

then  be  cancelled  and  Equation  B-22  reduces  to 


[(JTF)T/6(-2g  ACT  M2)T/6]T/2  = 


(t/6)  V  (ze  +  ze  t  1) 

3  Z-res  2,  .  w  3* 

ze(ze  -  i)(z2  -  ze) 


ze  =  o,i 


To  calculate  the  residue  at  the  ze  =  1  pole,  the  Ze  -  1  factor  is  deleted 
and  the  remaining  expression  is  evaluated  at  ze  =  1. 
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(B-24) 


(26  +  Zg  +  1)  =  (l2  +  1  +  1)2  =  9 

Z6(z2  -  zj)  12(Z2  -  l3)  Z2  '  1 

ze  =  1 


2 

Since  the  pole  at  zg  =  0  has  a  multiplicity  of  t\*o  (i.e.,  zg  =  0) ,  the 

2 

first  derivative  of  Equation  B-23  must  be  obtained  (with  the  zg  pole  deleted) 
before  calculating  the  residue.  That  is. 


(z|  +  26  +  1>2 

d  C(Z6  -  l)(z2  -  Z|)]  =Q 

dze 


(B-25) 


or 


d 

dzg 


(zj  +  ze  +  l)2 

C^2  -  4  -  *2  +  sl)]  Z6  =  o 


(B-26) 


Taking  the  first  derivative  in  Equation  B-26  with  respect  to  Z6  produces 
(Z6Z2  “  z|  -  Z2  +  z|)  [2(zg  +  Z6  +  l)(2zg  +  1)]  -  (z|  +  zg  +  l)2(z2  -  4zg  +  3z|) 


Adding  the  residues  from  Equations  B-24  and  B-28  gives  the  desired  low- 
rate  z-plane  transfer  function  for  item  1. 


(ITEM  1)  =[(JTF)t/6(-2_  ACT  M2)t/6]T//2  =  (T/18)  ( — - — T  -  — ) 

DE  Z2  -  1  z2 


=  (T/18) 


6z2  +  3 
z2(z2  -  1) 


(B-29) 


z2  =  esT/2 


(B-30) 


ITEM  2  (ATF)T/2(DTF)T/2C(JTF)T/6(-^  ACT  M2)T/6]T/2 


Substituting  the  appropriate  transfer  functions  from  Equation  B-5  and  the 
transfer  function  in  Equation  B-29  gives  ' 

<™2>  (T/18) 


(z2  +  1)2(6z2  +  3) 
z2(z2  -  l)3 


(B-31) 


where 

K  =  t3/(16)(18)  (B-32) 

ITEM  3  [1  -  (ATF)t/2(DTF)t/2C(JTF)t/6(-^  ACT  M2)t/6]T,/2]  1 

Substituting,  the  inverse  becomes 

ina  3>  -  ti  -  k  ia  i  .1)2(6Z2  -  —i 

Z2(z2  -  1)3 


-1 


— -v.  -i  ■  — r  ■ 


U 


m 


z2(z2  “  “  K  (z2  +  l)2(6z2  +  3)  1 

C  Z2(z2  -  1)3 


z2(z2  -  l)3 


z4  -  (3  +  6K)z3  +  (3  -  15K) z?  -  (1  +  12K)z2  -  3K 


(B-33) 


ITEM  4  Gj/2  =  [1  -  (ATF)t/2(DTF)t/2C(JTF)t/6(-2_  ACT  M2)t/6]T//2]_1(ATF)t/2 


DE 


Substituting  for  filter  (ATF)t/2  and  utilizing  Equation  B-33,  item  4  is 
given  by 


(ITEM  4)  = 


(T/4)z2(z2  -1 ) 2 ( z2  +  1) 


z^  -  (3  +  6K)z3  +  (3  -  15K) z?  -  (1  +  12K)z2  -  3K 


(B-34) 


ITEM  5  [(jTF)T/6(^  act  M2)T/6]T/2 


Using  the  same  procedures  as  for  item  1  results  in  the  following: 


(H  ACT  M2)T/6 


(1  -  e"57/2  l^T/6 


(1  -  Z^Hi^T/6  .  ill  V 

6  s3  z3  2(zg  -  l)3 


_  (zg  -  l)(z|  +  Zg  +  1)  (T/6)2z6(z6  +  1) 


z3 

Z6 


2(zg  -  l)3 
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I 


(zg  +  Zg  +  l)(zg  +  1) 
z|(zg  -  l)2 


(B-35) 


where 


KX  =  (T/6)2/2 


(B-36) 


[(JTF)t/6(~  ACT  M2)T/6]T/^2  = 


r  zj  +  *6  +  1  (zj  +  %  +  +  1)1  T^2 

3z|  1  z2(zg  -  l)2 


2  2  " 

r(z6  +  26  +  !)  (Z6  +  1), 

2C - zff^TTP - J 


(B-37) 


In  Equation  B-37, 


K2  =  Ki/3  =  (T/6)2/6 


(B-38) 


The  next  step  is  to  substitute  Equation  B-37  into  Equation  B-20  and  evaluate 
the  residues  associated  with  the  poles  at  zg  =  0  and  zg  =  1. 


(ITEM  5)  =  K2  7  residues 


(zg  +  zg  +  l)2(zg  +  1)  z2 


z^(zg  -  1)2 
6 


z2  -  z3 


(B-39) 


zg  *  0,1 
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Again  recognizing  that  zg  =  z2zg,  Equation  B-39  can  be  simplified  to 


(ITEM  5) 


residues 


(zg  +  zg  +  l)2(zg  +  1) 
z|(zg  -  1)2(z2  -  z|) 


zg  =  0,1 


The  residue  associated  with  the  pole  at  zg  =  0  is  given  by 


(zg  +  zg  +  l)2(zg  +  1) 
a  (zg  -  l)2(z2  -  zjb  . 

dzg 


zg  =  0 


5_ 

z2 


In  a  like  manner,  the  residue  for  the  Zg  =  1  pole  is  given  by 


d 

dzg 


(z|  +  zg  +  l)2(zg  +  1) 
c  - ] 


Zg  =  ! 


9(z2  +  5) 
(z2  -  l)2 


The  lew-rate  z-plane  transfer  function  form  item  5  is  obtained  by 
residues  in  Equations  B-41  and  B-42. 


(ITEM  5)  =  K2  +  9(Z2  +  5)]  =  K_  14z2  +  35z2  +  5 
Lz2  (z2  -  1)2J  *2  z2(z2  -  1)2 


ITEM  6  (ETF)T/2(ITF),r/2G^/2C(0TF)T/6(-^  ACT 


(B-40) 


(B-41) 


(B-42) 


adding  the 


(B-43) 


Substituting  Equations  B-5,  B-6,  B-34,  and  B-43  into  item  6  produces 


(ITEM  6)  =  (T/4) 


z2  +  1  _ (T/4)z2(z2  ~  1)2(z2  +  1) _ 

z2  ~  1  z*  -  (3  +  6K)z3  +  (3  -  15K)z2  -  (1  +  I2K)z2  -  3K 


14z2  +  35z2  +  5 
z2(z2  -  l)2 


-  k3 


_ (z2  +  1)2(14z2  +  35z2  +  5) _ 

(z2  -  1)  [z§  -  (3  +  6K)z^  +  (3  -  15K) z^  -  (1  +  12K)z2  -  3K] 


(B-44) 


Where 


K3  =  (T/4)2K2  =  (T/4)2(T/6)2/6  (B-45) 

ITEM  7  [1  +  (FTF)T/2(ITF)T/2GT/2[(jTF)T/6(^  ACT  M2)T/6]T/2]“1 


Utilizing  Equation  B-44,  the  inverse  in  item  7  is  given  by 


(nTM  7)  =  [1  +  k3  _ , _ _ a 

(z2  -  1)[ z£  -  (3  +  6K)z2  +  (3  -  15K) z2  -  (1  +  12K)z2  -  3K] 


-1 


4  3  2 

(z2  —  l)(z2  +  a3z2  +  a2z2  +  a^z2  +  a^) 

zf  +  b4z^  +  b3z|  +  fc^z2  +  biz2  + 


(B-46) 


where 


b4  *  -(4  +  6K  -  14K3) 
b3  =  (6  -  9K  +  63K3) 
b2  =  -(4  -  3K  -  89K3) 


-  • 
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a3  =  -(3  +  6K) 
a2  =  (3  -  15K) 

a1  *  -(1  +  12K) 


r— r- 


K 


M 


m 


'm 


\  *\ 
•  *.'t 

* 


ITEM  8 


<?/2= 


hi  -  (1  +  9K  +  45K3)  ao  =  -(3K) 

bo  -  (3K  +  5K3) 

K  =  T2/ (16) (18)  K3  =  (T/4)2(T/6)2/6 

Cl  +  (PTF)t/2(ITF)t/2gT/2[(jtf)t/6(-^  ACT  M2)T/6]T/2f1 


*  (ftp)t/2^/2 


Before  forming  G^/2,  the  transform  N^1/2  must  be  calculated.  Fran  Equation 
B-9  we  have 


dj/2  =  (■ 


.  -ST  T/2  ,  T/2 

1  -  (i  -  *ix±> 


z2  ~  1  z2  z2  +  1 

z|  z2  ~  1  ~  z2 


(B-47) 


Combining  Equations  B-46  and  B-47  with  the  (FTF)T/2  filter  in  Equation  B-5 
gives 


2#_4  .  _  3  .  _  _2 


(ITEM  8)  =  d£/ 2  =  (T/4)(*2  +  1)  <z2  +  a3z2  *  *2Z2  +  alz2  +  ap> 

“  z3(z2  +  ^4*2  +  ^3z2  +  ^2z2  +  ^lz2  +  ^0^ 


item  9  ^  =  [(||  PCT  M2)T/W2f  GC 


The  s-plane  to  z-plane  transform  in  item  9  is  given  by 
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{1  “  z51)(^r) 

z2  -  1  (T/2)2z2(z2  +  1) 

z2  2(z2  -  l)3 


(T2/8)(z2  +  1) 
(z2  -  l)2 


(B-48) 


=  (t2/8)(z2  +  1)  (T/4)z2(z2  -  1)2(z2  +  1) 

A?  ^  <*2  -  D2  (z^  +  a3z^  +  a2z2  +  aiz2  +  bq) 


(T/4)(z2  +  1)2(z2  +  a2z2  +  a2z2  +  a^z2  +  a q)  i  1 
z2(z3  +  b4z^  +  t^z3  +  l^z2  ^  b1z2  +  b^)  * 


c _ (TVl28)(Z2  -f  l)4 _ 1 

(z3  +  b4z^  +  tgz3  +  l^z2  +  b^z2  +  fc^) 


(&-49) 


Hie  evaluation  of  the  high- to- lew  rate  transform  conversion  in  Equation 
B-49  is  tedious.  The  results  of  applying  option  149  of  TOTAL  is  shewn  belcw. 


aS  =  (-.8024E-13)z5  +  (.2441E-0DZ4  +  { .1709E-06)z3  +  ( .1709E-06)z2 

ASc  _ -f  ( .2441E-07)z  -  ( .9968E-15) _ 

(1.000)z5  -  (4.000)z4  +  (6.000) z3  -  (4.000)z2  +  (1.000)z 


-  (.2667E-13) 


(B-50) 


Equation  B-50  is  the  desired  closed-loop  z-plane  transfer  function  defined 
for  z  =  eP?. 


.  ‘  ■  %  ••  ■*  /.  /.  *.  ••  A  A 
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APPENDIX  C 


MECHANIZATION  OF  MULTI-RATE  MODEL  EQUATIONS  USING  TOTAL 

INTRODUCTION 

An  automated  technique  for  evaluating  multi-rate  model  equations  is 
presented.  This  technique  enploys  the  transfer  function  storage  and  macro 
options  in  the  linear  analysis  computer  program  TOTAL  (Reference  13) .  The 
closed-loop  multi-rate  model  equation  in  Section  V  (Equation  103)  is  used  to 
illustrate  this  technique.  Equation  103  is  repeated  below  as  Equation  C-l. 

4  -  C(|g  KT  l+,)T/2GT/2GT/2]T  (C-l) 

AZq 

In  Equation  C-l,  the  following  definitions  apply  (see  Equations  104-106) 

<3T/2  »[l  -  (ATF)T/2(ETF)T/2C(JTF)'r/6(^  ACT  M2)T/6]T/2]~1  (ATF)t/2  (C-2) 

G|T/2  *  [!  +  (rTF)T/2(lTF)T/2GT/2C(JTF)T/6(-^  ACT  M2)T/6]T/2]  1 

*  (FTF)T/2^/2  (C— 3) 

GT  =  (KTF)T  (C—4) 

To  apply  this  automated  technique,  the  s-plane  and  z-plane  transfer 
functions  that  appear  in  Equations  C-l  through  C-4  sure  stored  in  TOTAL.  A 
user  generated  macro  is  then  executed  that  accesses  these  transfer  functions 
and  performs  a  sequence  of  operations  to  solve  the  multi-rate  equation.  The 
transfer  function  storage  and  macro  capabilities  of  TOTAL  are  described  in 
the  following  subsections. 


TRANSFER  FUNCTION  STORAGE 


TOTAL  has  24  auxiliary  transfer  function  storage  locations  designated  by 
ATF,  BTF,  ...  ,  ZTF  (GTF  and  HTF  serve  other  functions  within  the  program 
and  are  excluded  fran  this  list) .  This  specific  transfer  function  notation 
is  used  in  the  multi- rate  model  equations.  However,  in  Equations  C-l  through 
C-4,  storage  locations  must  also  be  assigned  to  N^/2,  ,  ACT,  q/DE,  and  AZ/DE. 

In  Equation  C-3,  M£/2  is  defined  by 

T/2 

1^/2  .  (1  Y~")  (05) 

To  evaluate  Equation  C-5,  the  z-plane  factor  is  pulled  out  and  the  T/2 
transform  is  obtained  for  the  remaining  s-plane  factor.  That  is, 

rtp  =  <1  -  z52)(i>T/2 


2  . 

z2  -  1  z2 
z2  z2  "  1 


Z2  +  1 
z2 


where 


Z2  =  esT/2 


(C-6) 


(0-7) 


For  the  macro  calculations,  we  let 

f^/2  =  (LTF)T/2  =  (c_8) 


The  M2  data  hold  appears  in  the  following  three  transforms: 


(|§  ACT  M2)T/2 
ACT  M2)T/6 
#  ACT  M2)T/6 


(C-9) 


(C-10) 

(C-ll) 

In  Equation  C-9,  it  is  not  necessary  to  explicitly  store  M2  in  TOTAL  to 
calculate  the  T/2  9-plane  to  z-plane  transform.  Since  M2  is  referenced  to  a 
T/2  sampling  interval  (M2  =  1  -  e_sT/2/s) ,  option  141  in  TOTAL  will  insert 
M2  before  taking  the  transform.  The  user  simply  supplies  the  Az/DE  and  ACT 
s-plane  transfer  functions. 

Hie  above  simplification  for  the  M2  data  hold  is  not  possible  in 
Equations  C-10  and  C-ll  since  the  sampling  interval  is  T/6.  It  will  be 
necessary  to  separate  M2  into  z-plane  and  s-plane  factors  and  place  each  in 
a  separate  transfer  function  storage  location.  Substituting  for  M2, 
Equations  C-10  and  C-ll  become 

t/s 

&  **  m2>t/6  -  ^ 1  ~  r/2> 


3 

_  *6 


-  1 


z3 

Z6 


T/6 


(C-12) 


<w  «2>t/6  -  #  1  -  r/2 


T/6 


3  ,  T/6 


(C-13) 


where 


(C-14) 


zg  =  e^/6 

For  the  macro,  the  z-plane  and  s-plane  factors  in  Equations  C-12  and  C-13 
will  be  stored  as 

3  _ 

(MrF)T/6  *  — = —  (C-15) 

Z6 

NTF  =  i  (C-16) 

S 

The  remaining  transfer  function  storage  allocations  are  given  below. 


BTF  =  ACT 

(C— 17) 

(71’F  =  — ' 

r  DE 

(C-18) 

.  Az 

(C— 19) 

YTF  *  y 

(C-20) 

ZTF  = 

(C-21) 

The  unity  transfer  functions  YTF  and  ZTF  are  needed  in  the  macro  to  perform 
the  inverse  calculations . 

Utilizing  the  above  storage  allocations  for  M^/2,  tfyi  ACT,  q/DE,  and 
Az/DE,  Equations  C-l  through  C-4  can  be  expressed  as 

T 

*  C(ETF  EfTF  M2)T/2GT/2GT/2]T  gT  (C-22) 

ASc  ABC 


G?/2  «  [i  _  (ATF)T/2(OTF)T/2[(JTF)T/6(MrF)T/6(CTF  BTF  otf)T/6]T/2]  1 


C^/2  =  [i  +  (PTF)T/2(lTF)T/2Gj/2[(i;rrF)T/6(MTF)T/6(ETF  grp  otF)T/6]T/2-j  -1 


*  (FTF)t/2(LTF)t/2  (C-24) 

<JT  =  (KTF)T  (C-25) 

As  pointed  oat,  it  is  not  necessary  to  store  the  M2  transfer  function  that 
appears  in  Equation  022.  It  will  be  inserted  by  option  141  in  TOTAL 
before  the  T/2  s-plane  to  Z2~plane  transformation  is  calculated. 

MACRO  INSTRUCTION  SET 

The  macro  capability  of  the  TOTAL  computer  program  can  be  used  to  per¬ 
form  the  series  of  operations  in  Equations  022  through  025 .  A  macro 
operates  as  an  executive  instruction  set  far  TOTAL.  Each  macro  contains  up 
to  98  sequential  instructions.  Typical  instructions  include  the  following: 
(1)  copy  a  transfer  function  from  one  storage  location  to  another,  (2) 
multiply  or  add  two  transfer  functions,  (3)  set  an  internal  variable  to  a 
numerical  value,  (4)  set  a  logic  switch  to  true  or  false,  (5)  perform  a  13- 
plane  to  z-plane  transformation,  and  (6)  calculate  a  low-rate  z-plane 
transfer  function  from  a  given  high-rate  z-plane  transfer  function.  These 
instructions  represent  a  small  sanple  of  the  large  number  of  options 
available  in  TOTAL  (there  are  over  100  active  options) . 

The  first  step  in  evaluating  the  multi-rate  equations  is  to  sequentially 
list  the  necessary  individual  operations.  Each  operation  is  then  related  to 
a  specific  series  of  catmands  and  options  in  TOTAL.  Finally,  by  invoking 
the  special  CREATE  ocrrmand  in  TOTAL,  this  sequential  list  of  carmands  and 


options  is  converted  to  a  coded  macro  instruction  set. 

The  AKEY,  BKEY,  and  CXEY  macros  for  Equations  0-22  through  C-25  are 
given  in  Tables  C-l ,  C-2 ,  and  C-3.  Column  one  in  these  tables  sequentially 
lists  the  individual  operations,  and  colurtn  two  lists  the  corresponding 
TOTAL  commands  and  options .  The  sanpling  intervals  used  in  these  macros  are 
T  =  1/40,  T/2  =  1/80,  and  T/6  =  1/240.  Table  C-4  outlines  the  specific 
TOTAL  ccnmands  and  options  utilized  in  the  macros. 

As  specified  by  the  user,  TOTAL  stores  the  coded  macro  instructions  in 
local  files  AKEY,  BKEY,  or  CXEY.  These  local  files  can  be  cataloged  as  per¬ 
manent  files  for  future  use.  For  Equations  C-22  through  C-25,  it  is 
necessary  to  divide  the  instructions  between  the  three  macros.  Each  macro 
is  limited  to  98  instructions,  and  232  instructions  are  required  to  evaluate 
these  multi-rate  equations.  The  AKEY  macro  contains  93  instructions,  BKEY 
85  instructions,  and  CKEY  54  instructions. 

It  is  the  users  responsibility  to  keep  track  of  the  instruction  count. 

In  each  macro,  a  complete  copy  operation  takes  three  instructions  (e.g. , 
COPY,  NTF,  GTF) ;  setting  a  variable  to  a  numerical  values  also  takes  three 
instructions  (e.g.,  TSAMP  =  .004167);  and  an  option  number  requires  one 
instruction  (e.g.,  21,  24,  141,  or  149). 

To  execute  the  AKEY  macro  in  TOTAL,  the  ccnmand  AKEY  is  given  while  in 
the  option  mode.  That  is, 

OPTION  >  AKEY 

TOTAL  will  then  go  to  the  local  file  AKEY  and  read  and  execute  each  instruc¬ 
tion  in  a  sequential  fashion.  From  Table  C-l,  the  first  three  instructions 
will  copy  the  transfer  function  stared  in  NTF  into  the  GTF  storage  location. 


Following  this,  TOTAL  will  copy  the  contents  of  BTF  into  ffTF.  Option  21 
will  next  multiply  GTF  and  HTF  and  store  the  resulting  transfer  function  in 
OLTF.  The  last  instruction  in  AKEY  is  to  call  BKEY.  When  TOTAL  encounters 
this  ccumand,  it  iimediately  begins  to  read  the  instructions  in  the  local 
file  BKEY.  In  a  similar  manner,  the  last  instruction  in  the  BKEY  macro  is 
to  call  the  CXEY  macro. 


— - - - — - - - - — 

CN 

• 

TABLE  C-l . 

AKEY  MACRO  - 

OPERATIC*! 

command/ option 

BTF  OTF 

COPY,  OTF,  GTF 

COPY,  BIF,  BTF 

• 

21 

CTF  BTF  OTF 

COPY,  OLTF,  GTF 

COPY,  CTF,  BTF 

J 

21 

• 

GT/6  -  (Grp  grp  ntF)T/6 

COPY,  OLTF,  CLTF 

TSAMP=. 004167 

H0LD=4 

DELTAfO 

141 

• 

(MTF)t/6gT/6 

COPY,  CLTF,  GTF 

COPY,  MTF,  BIF 

21 

(  JTF)t/ 6  (MTF)  T/6gT/6 

COPY,  OLTF,  GTF 

• 

COPY,  JTF,  BTF 

21 

<5/2  =  [(JTDT/SCMTFjT/eGT/e]^ 

COPY,  OLTF,  CLTF 

HISAMP=1 .0 

-  '  .. 

LOSAMP=3 .0 

149 

• 

(OT F)T/2gT/2 

COPY,  CLTF,  GTF 

•  •  .  ■ 

COPY,  OTF,  BTF 

21 

/ 

• 

gT/2  =  (OTF)t/2(DTF)t/2gT/2 

QDPY,  OLTF,  GTF 

COPY,  ATF,  BTF 

21 

-<5/2 

COPY,  OLTF,  GTF 

o 

COPY,  ZTF,  BTF 

21 

• 

I  -  <5/2 

COPY,  OLTF,  GTF 

J 

COPY,  YTF,  BTF 

24 

• 

(I  -  (5/2)-1 

COPY,  CLNPOLY,  GDPOLY 

COPY,  CLDPOLY,  GNPOLY 

<^/2  =  (I  -  <^/2)_1(ATF)'r/2 

COPY,  ATF,  BTF 

21 

STORE  <5/2  IN  OTF 

A 

_  _• 

COPY,  OLTF,  CTF 

CALL  EKEY 

87 

EKEY 

.  • 

«-*  *-■’  *-•  fcliJ. ■-*  *-**»  *  ^  *  '  ii 

TABLE  02.  BKEY  -  G£/2 


OPERATION 

command/ option 

BTF  NTF 

COPY,  OTF,  GTF 

COPY,  BTF,  OTF 

21 

ETF  ETF  NTF 

OOPY,  OLTF,  GTF 

COPY,  ETF,  HTF 

21 

G?/6  =  (etf  grp  NTF)t/6 

COPY,  OLTF,  CLTF 

141 

(MTF)t/6gT/6 

COPY,  CLTF,  GTF 

OOPY,  MTF,  HTF 

21 

(  JTF)t/6  (MTF)T/6GT/6 

OOPY,  OLTF,  GTF 

COPY,  JTF,  OTF 

21 

GT/2  =  C(JTF)t/6(MTF)t/6gT/6]T/2 

COPY,  OLTF,  CLTF 

149 

GT/2QT/2 

A  2 

OOPY,  CLTF,  GTF 

OOPY,  OTF,  OTF 

21 

(lTF)T/2GT/2GT/2 

COPY,  CLTF,  GTF 

COPY,  ITF,  OTF 

21 

Gff/2  =  (FTF)T/2(ITF)T/2GT/2GT/2 

n  £ 

COPY,  OLTF,  GTF 

COPY,  FTF,  "OTF 

21 

I+<?/2 

OOPY,  OLTF,  GTF 

OOPY,  YTF,  OTF 

24 

(I  +  dT/2)-1 

OOPY,  CLNPOLY,  GDPOLY 
OOPY,  CLDPOLY,  GNPOLY 

(I  +  <^/2)“1(PTF)T/2 

OOPY,  FTF,  OTF 

21 

G^/2  =  (i  +  GT/2)"1(FTF)T/2(LTF)T/2 

OOPY,  OLTF,  GTF 

OOPY,  LTF,  OTF 

21 

STORE  dj/2  in  PTF 

D 

OOPY,  OLTF,  PTF 

CALL  CKEY 


CXEY 


2/2 


P 


AD-A132  844  MULTI-RATE  DIGITAL  AUTOPILOT  MODELING  HETHODOLOGV(U) 

ARMAMENT  DIV  CAFSC)  EGLIN  AFB  FL  D  G  DIDALEUSKV  MAR  83 
AD-TR-83-64  AFATL-TR-83-S7 

UNCLASSIFIED  F/G  1/2 


NL 


OPERATION 


TABLE  0-3.  CXEY  MACRO  -  aI/AZc 

CXMttND/OPTION 


ETF  BTF 


dj/2  *  (ETF  EfTF  M2)T/2 


g3’/2GT/2 

1  A 


g3’/2GT/2gT/2 

1  A  d 


(GT/2GT/2GT/2)T 


^  -  «^/2gT/2gT/2)TST 

Azc 


STORE  Az/AZc  IN  QTF 
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GOPY,  BIT,  GTF 
COPY,  ETF,  HIT 
21 

COPY,  CLTF,  CLTF 

TSAMP=.0125 

HQLD=0 

DELTA=0 

141 

COPY,  CLTF,  GTF 
COPY,  OTF,  HIT 
21 

COPY,  CLTF,  GTF 
GOPY,  PTF,  HTF 
21 

COPY,  CLTF,  CLTF 
HISAMP=1.0 
LOSAMP^2.0 
149 

COPY,  CLTF,  GTF 
GOPY,  KTF,  HIT 
21 

COPY,  CLTF,  QTF 


•  .*  A  . 


.  /»  / 


TABLE  C-4.  TOTAL  OOWffiNDS  AND  OPTIONS 


QOWANDS/OPTIONS 
COPY,  OTF,  GTF 

21 

24 

141 

TSAMP 

HOLD 

DELTA 

149 

HISAMP 


DESCRIPTION 


Copies  transfer  function  stored  in 
OTF  into  storage  location  GTF. 
Leaves  contents  of  OTF  unaltered. 

Option  21. 

GLTF  *  (jTF  *  OTF 
Multiplies  transfer  functions 
stored  in  GTF  and  OTF  and  stores 
product  in  OLTF. 

Option  24. 

CLTF  =  GTF  +  OTF 

Adds  transfer  functions  stored  in 
GTF  and  OTF  and  stores  sun  in  CLTF. 

Option  141. 

Converts  s-plane  transfer  function 
stared  in  CLTF  to  z-plane  transfer 
function  and  stares  results  in  CLTF 

Sampling  interval  (sec). 

Hold  option:  (0)  zero-order,  (1) 
first-order,  (2)  second-order,  (3) 
slewer,  or  (4)  none. 

Transform  option:  (0)  standard 
transform,  (+)  advanced  transform, 
or  (-)  delayed  transform. 

Option  149. 

Converts  high-rate  z-plane  transfer 
function  stored  in  CLTF  to  lew-rate 
z-plane  transfer  function.  Stares 
results  in  CLTF. 

High  rate  sampling  interval. 

Lew  rate  sampling  interval  (Only 
ratio  of  sampling  intervals  needed 
for  HISAMP  and  LOSAMP.) 


LOSAMP 


INITIAL  DISTRIBUTION 


DTIC-DDAC  2 

AUL/LSE  1 

FTD/SDNF  1 

HQ  USAFE/INAT  1 

AFWAL  (FIES/CDIC)  1 

AFWAL  (FIGL)  2 

COMDR  NAVAIRSYSCOMD  (CODE  AIR-5401E)  10 
AD/ENMA  IQ 

AD/YMEB  3 

AD/XRCD  1 

AFATL/DLMA  2 

LEAR  SIEGLER,  INC.  1 

GEN  ELECTRIC  COMPANY  1 

BOEING  AEROSPACE  COMPANY  1 

NORTHROP  CORPORATION  1 

WRIGHT  STATE  UNIVERSITY  1 

AFALT/DLODL  2  * 
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